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FOREWORD 


In the past fifteen years the wave equation has 
come to have a steadily increasing importance in theoretical 
physics. In addition to the increase in practical impor- 
tance of the classical fields of acoustics and electro- 
magnetic theory, the new and fundamental field of wave 
mechanics has come into being. The developments in this 
field have shown that a basic understanding of the structure 
of matter and of particle dynamics can only be obtained 
through a study of a wave equation. Thus, both old fields 
of application of wave theory and basically important new 
fields have indicated an expanding need for an improvement 
in the techniques for the calculation of solutions of the 
wave equation. 

For the solution of many general problems in- 
volving wave motion, general solutions of the wave equation 
can be used. For instance, in many cases a sum, or inte- 
gral, of plane waves may be used to compute general behav- 
ior; although in many such computations serious questions 
of convergence of the series may arise. Recent calcula- 
tions in quantum electrodynamics are typical examples of 
the plane wave summation technique, and also of the diffi- 
culties this technique encounters. There is increasing 
need, however, for the solution of specific and detailed 
problems rather than general cases; and for this work 
representations of the wave function must be found which 
are better- and more rapidly convergent than the plane wave 
aggregates. 

This need for better representations is partic- 
ularly apparent with wave problems involving boundary 
conditions on finite boundaries. In such cases, usually 
the only practicable method of numerical solution involves 
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thG separation of the wave equation in the particular co- 
ordinate system which corresponds to the boundary. Thus 
the partial differential equation is reduced to a sequence 
of ordinary differential equations, one for each co-ordinate. 
This method is exceedingly powerful in allowing specific 
boundary value problems to be solved in detail j but it is 
severely limited in its scope, for the wave equation is 
separable in only a few Euclidian co-ordinate systems. The 
work of Robertson^ and of Eisenhart^ has shown that only 
eleven different co-ordinate systems allow separation of 
the scalar wave equation in three dimensions. Some further 
work has indicated that the limitations for the vector wave 
equation in this respect are even more binding, for it ap- 
pears likely that only five different tri-dimensional co- 
ordinate systems allow vector solutions of the wave equation 
which can be separated in such a manner as to satisfy 
boundary conditions on the co-ordinate surfaces. 

This limitation of the applicability of the separa- 
tion technique is a serious one and it appears urgent that 
further fundamental research be done in the development of 
other methods of solution. Nevertheless, from the point of 
view of the application of wave theory to physics and to 
engineering, the possible usefulness of the separation tech- 
nique is far from having been exhausted. In order to apply 
this technique to a particular problem the solutions of the 
ordinary differential equations arising from the separation 
must be obtained, their mathematical properties must be 
thoroughly explored, and tables of values of the solutions 
must be prepared. Only after these three steps have been 
taken can numerical solutions to specific problems be com- 
puted without great difficulty. This state of affairs has 
been reached at present for only three of the eleven co- 
ordinate systems: rectangular, circular cylinder, and 


vi 



spherical co-ordinates- 

The ordinary differential equations arising from 
the separation of the wave equation, in the eleven cases 
mentioned, all turn out to be second order and linear, with 
the singular points of the equation determined primarily by 
the geometrical properties of the particular co-ordinate 
system. For instance, the equations arising from the separa- 
tion into rectangular co-ordinates have a single irregular 
singular point at infinity; whereas some of the equations 
for the circular cylinder and spherical co-ordinates have 
one regular singular point at zero and one irregular point 
at infinity, and other equations have three regular singu- 
lar points. The mathematical properties of the solutions 
of equations with two and three regular singular points, one 
irregular, or one regular and one irregular point, have been 
quite thoroughly studied in the past century. Tables of 
values of some of the solutions (trigonometric functions, 
Bessel functions, spherical harmonics, etc.) have been com- 
puted and are available in various states of completeness. 
Thus the practical problem of solutions for plane, cylin- 
drical, or spherical waves has been made as easy as possible 
for the engineer and applied physicist. The solutions for 
the other eight separable co-ordinate systems, however, are 

in no such satisfactory state. 

In order of practical importance, the next three 

co-ordinate systems to be made available for application 
should be the elliptic cylinder, the prolate spheroidal, and 
the oblate spheroidal co-ordinates. Solutions of problems 
involving the radiation and scattering of waves from strips 
of material, from wires of finite length, and from discs of 
material, ail require the Knowledge of the mathematical 
properties and the numerical values of solutions of the wave 
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equation for these three co-ordinate systems. The solutions 
are likewise required for the study of the diffraction of 
waves through slits and circular openings, the absorption of 
sound by strips or by circular patches of material, and the 
behavior of electrons in diatomic molecules, to mention a 
few cases. 

The ordinary differential equations arising from 
the separation of the wave equation in these three co- 
ordinate systems have two regular and one irregular singular 
points. The mathematical properties of the solutions of 
equations of this sort have been studied for some time. 
Mathieu^ began the study of their application, with his 
classical study of the vibrations of an elliptical membrane, 
since that time Ince^, Humbert 5 , Strutt 6 and many others 
have made contributions. A consolidation of the knowledge 
of these solutions has not been completed as yet, however. 

In particular, a standard form of solution and a standard 
notation for the solutions has not yet been agreed upon. 
Since such agreement in other cases has usually accompanied 
or followed the publication of reasonably complete tables 
of solutions, it is hoped that the present volume will help 

to make agreement possible for elliptic cylinder and 

7 

spheroidal functions. 

The present volume defines certain standard forms 
of the solution which turn out to be of use in practical 
problems, and it displays a collection of formulas giving 
the important mathematical properties of these functions. 
Many of these properties have been discovered and studied 
earlier. They are given here in terms of the standard form 
of solution in order to make the outline complete. 

In addition to the mathematical properties of the 
solutions, there is also a set of tables from which values 
of the solutions can be obtained for the more interesting 
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ranges of the variables. These tables contain the series 
coefficients together with some of the separation - constants 
and allied coefficients, and cover the most difficult part 
of the whole procedure of obtaining numerical solutions. 
Although it is expected that eventually a larger volume will 
be published which will give actual numerical values of the 
solutions themselves, nevertheless, these numerical solu- 
tions can be obtained relatively quickly from the present 
tables of coefficients. 

It is to be hoped that the present volume will 
begin to make it possible for the applied physicist and 
engineer to handle wave problems in elliptic-cylindrical 
and spheroidal co-ordinates with approximately the same 
degree of facility as has been possible previously for rec- 
tangular, circular cylindrical, and spherical co-ordinates. 

The present tables are the result of a number of 
years of effort on the part of the authors. Many persons 
at the Massachusetts Institute of Technology have contrib- 
uted to the work and the complete list of acknowledgments 
would be a long one. Miss Pearl Rubenstein helped one of 
the authors publish a preliminary table of Mathieu functions 
several years ago. The present tables are to some extent 
based on these earlier ones. The N. Y. A. has contributed 
toward the payment of a number of students in order to help 
in some preliminary calculations. The Department of Physics 
at the Massachusetts Institute of Technology has contributed 
computing machines and other equipment, and the Differential 
Analyzer staff has been of considerable aid in the work. 

The writers wish to express their appreciation to Professor 
S. A. Caldwell, in charge of the Massachusetts Institute of 
Technology Computing Center, for his help in this respect 
and to Mr. J. R. Killian for his valuable assistance in 
making possible the final calculations and the publication 



of these tables. I he authors also wish to express their 
appreciation of the continuing interest of President 
K. T. Compton in their lengthy task. 


Philip M. Morse 


Cambridge, Mass. 
October 18, 1941® 
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7„ It is desirable to call the reader's attention to the 
fact that in this volume the notation for the elliptic cyl- 
inder functions, and both the notation and the method of 
normalization for the spheroidal functions are new, and 
therefore different from those appearing in the previous 
papers published by any of the authors. 


For example, the function Se^(c,z) in the paper by 
Stratton, Proc. Nat. Acad. Sci., 21, 51, 1935, is now de- 
noted by v( X jJ(c,z), and the function ( c > z ) b y U^*(c,z), 

with similar ' changes for the solutions of the second kind. 

In the case of the elliptic cylinder functions, Se _ / ^ i ( c > z ) 

is now represented by Se^^Cjz), R e _^A ( c > z ) b y 


and He ^ \ (c,z) by Ne. (c,z), with similar changes for the 
odd functions. 


On the other hand, the formulas in spheroidal coordinates 
given in the paper by Morse, Proc. Nat. Acad. Sci., 21, 56, 
1935, are affected by the change in normalization. Thus 
equation (14) for the plane wave addition formula becomes; 

e ikX _ 2^Tcos [m(<p-a)] X d£ Xkt2ffiXL"j y 

mn 2 m mJ Lk=0,l J 

(° » cos “)'t r ^' ) (*} » 

where (c,cos0)R^ (c,coshn); 

the direction of propagation with respect to the z axis is 
defined by the angles co and a; and 

X = ^[cosco cosQ coshjA + since sin© sinhp cos(cp-a)] ; 

k = 2a g c = tfd, 

\ X 

However, in the elliptic cylinder coordinates, the 
normalization has remained unchanged. Thus the equations 
given in Morse's paper are still valid, though there is a 
slight change in notation. For example. 



e iM = T§I I i n 


n 


^ e ^ 1} (c,cosu)^ 1) (k;cp,ii) 

L n 


+ So^ 1 ^ (cjcosujx^ 1 -* (k;{p,|i)l 

N n J 


v/here i}^ 1 ^ (k;cp,|j,) = Se^ 1 ^ (c,cosc^)Je n (c, cosily) 5 

= So^ 1 ^ (c,cos(p)Jo n (c,cosh[i); 

u is the angle with respect to the x axis of the direction 
of propagation of the plane wave; and 

Y = — (cosu coscp coshp + sinu sincp sinhp) . 
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ERRATA 


p. 3 Equation (14): for "VxVx s" read "/uVx Vx s." 

p. 18 Top of page: for "T^ 11 read n rJ +r " 

p. 49 Equation (310): for "n \[E " read "t " 

r 7T f 2 

p. 54)57 First and second equations below "2. Radial Functions”': 

for "m" in the arguments of the cos and the sin read "J." 

p. 56 The equation for Aj ) L = even: for read "Dq." 

p. 57 Bottom of page: for "z = cos ^ 11 read u z = cosh^." 

p. 58 All E's and L's should have subscripts "jL.” 

p. 62 Third equation from the top: for "z = t? = £" read 

"z = Tj or £ . " 

p. 63 Top equation for n dP^ +n n read "dP^^ . ” 

?• 72 Third line: omit the word "negative." 

Thirteenth line: for "negative coefficients" read 
"coefficients with negative subscripts." 
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I. Introduction 

1.1 Scalar Wave Equation. It is known that the wave equation 1 

1 The most complete account of earlier researches on the subject matter of 
this paper is contained in a monograph by M. J. O. Strutt, “Lamesche, 
Mathieusche, und verwandte Funktionen in Physik und Technik,” Springer, 
1932. An extensive bibliography is included. 
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V 2 W + k 2 W = 0 (1) 

may be separated in spheroidal coordinates £, y, <p obtained from the 
cylindrical coordinates r, <p; z through the transformation 

r = fV(? + 1)(1 — y 2 )> <P = <P, z = f&, ( 2 ) 

the upper or lower sign being chosen according as the coordinates are of 
prolate or oblate type. The factors of the resulting Lam6 product are 


then 

m m. 

W = (£ 2 + 1) T um 1 - V(y)<H<p) (3) 

and satisfy the ordinary equations 

$"(*>) + = 0 , ( 4 ) 

(t* + 1 )U" + 2(m + l)£t/' - (b - fh 2 i?)U = 0, (5) 

(1 - y)V" - 2 (m + l)yV' + (b + AY) 7 = 0, (6) 

where b is a separation constant to be determined from conditions of 
finiteness and periodicity over the surface of a spheroid. Both (5) and 
(6) can be reduced to the standard form 

(1 - z 2 )U " - 2(o + 1 )zU’ + (6 - cz) U = 0. (7) 

The coordinates of the elliptic cylinder are expressed in terms of rec- 
tangular coordinates by the transformation 

V = fVV? - 1)(1 - -I 2 ), z = z (8) 

In this system the scalar wave equation (1) separates into the factors 

W = U(£)V(y)Z(z) (9) 

satisfying the ordinary equations 

Z" + y 2 Z = 0, (10) 

(f - 1)1/" + W - [b - f\k 2 - i*)t]U = 0 (11) 

(1 - v 2 )V" - -nV' + [6 - - t *)v]V = 0. (12) 


Equations (11) and (12) are obviously special cases of (7) for a = — J. 
Their solutions are known generally as Mathieu functions. 

To visualize the geometrical significance of certain limiting cases the 
following will prove helpful. The eccentricity of an ellipse £ = constant 
is e = l/£. Let l be the length of a semi-major axis. Then the coordi- 
nate £ of an ellipse is 
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£ = Z// = 1/e, (13) 

thus appearing to be a pure number, the ratio of semi-major axis to 
semi-focal distance, or reciprocal of the eccentricity. The limiting case 
of a circle occurs when / — > 0, £ — * ► °° such that the product £/ remains 
finite and equal to the limiting radius 1. Identical relations hold for the 
spheroids obtained by rotating an ellipse about a principal axis. 

1.2 Vector Wave Equation. In the theory of elasticity, as well as 
in electromagnetic theory, one must deal frequently with vector rather 
than scalar wave functions. Thus the equation governing the propaga- 
tion of elastic waves in a homogeneous, isotropic medium is 

d z s 

VxVss - (A + 2 m)VV-s + p -r = 0, (14) 

av 

where s is the relative displacement of a point in the medium, p is the 
density and A and p are elastic constants. Likewise an electromagnetic 
field whose angular frequency is co can be derived from the Hertz vectors 
II and IX* according to the rules 2 

E = k 2 n -f V* — iupV * n* (15) 

H = (ia>€ + <r)V X II + tfll* + V^*, (16) 

k 2 = peJ — ipcru, (17) 

where p, e, a are respectively permeability, dielectric constant, and 

conductivity in rationalized m.k.s. units. The electric and magnetic 
Hertzian vectors II and II* are independent and satisfy 

-Vx VxII + A; 2 XI + Vyp = 0. (18) 

The scalars i p and \p* are likewise independent. Usually they are placed 
equal to the divergence of the corresponding Hertz vector, but any solu- 
tion of a homogeneous scalar wave equation may be employed. 

Solutions of vector equations such as (14) and (18) can be constructed 
from solutions of the scalar -wave equation. 3 Thus any rectangular 
component of a Hertzian vector satisfies (1). The differential operators 
may be expressed in terms of curvilinear coordinates and upon carrying 
out the operations indicated in (15) and (16) an electromagnetic field 
is obtained in a curvilinear system of coordinates. Such fields must 

2 J. A. Stratton, “Electromagnetic Theory,” pp. 28-32, McGraw-Hill, 1941. 

3 J. A. Stratton, loc. cit., pp. 349-354, 392-399, 414-420. 
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eventually be superposed to satisfy boundary conditions on a prescribed 
surface. In place of a rectangular component it may be shown that the 
spherical function Ur/R is also a solution of the scalar wave equation, 
where II a is a radial component. The field derived from a single com- 
ponent n« has no radial component of H; that derived from 11* has no 
radial component of E. Unfortunately a vector field satisfying boun- 
dary conditions on a spheroid cannot be obtained in the same way from 
Hertz vectors which are everywhere normal to the surface of the 
spheroid. Such fields must be constructed by superposition of partial 
fields derived in the manner just described. 

1.3 Gegenbauer Functions. The standard equation (7) is charac- 
terized by regular singularities at z — ±1 and an irregular point at 
infinity. The object of this investigation is to obtain independent solu- 
tions in the neighborhood of each singularity reduced to a form suitable 
for numerical calculations and to discuss their analytic continuation 
over the entire complex plane of z. 

From (13) we note that if / — » 0, £ — » «> in such a way that /£ — > R, 
the spheroid reduces to a sphere of radius R. In this limit (7) is re- 
placed by 

iU" + 2(o + 1 )zU' + (eV - b)U = 0, (19) 

an equation satisfied by half-order Bessel functions. For this reason 
the spheroidal functions will be represented in the neighborhood of the 
essential singularity at infinity by expansions in Bessel functions. 

If, on the other hand, we place c = 0, (7) reduces to 

(1 - z 2 )U" - 2 (a + 1 )zU' + W = 0. (20) 

This is essentially the Gegenbauer equation, a hypergeometric equation 
with three regular singularities. 4 Since the spheroidal functions will be 
expressed in terms of the solutions of (20), we shall give a brief account 
of their properties. 

The functions T a n (z) and ( z ) are defined for all values of a, n, and 

z by 

n« = (*» - i )-" r, +a (z), (2i) 

r.w = («* - ir n Qi+M, (22) 

where P a n+a (z) and Q a n + a (z ) are the associated Legendre functions of the 

4 Whitaker and Watson, “Modern Analysis,” p. 329, 4th ed., Cambridge Uni- 
versity Press. 
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first and second kinds. Definitions of the Legendre functions differ in 
the literature with regard to phase, or the location of the cut in the 
2 -plane. Accordingly, w§ shall adhere to the definitions given by 
Hobson. 6 


K(z) = 


e nrt (n -j- a ) ! ^2 __ ^<V 2 
4t sin tot 2 n n\ 


X / « 2 

J C 


l) n (t - z)- n - a ~Ut (23) 


Q a n(z) 


i e n1ti (n 4~ a) 1 ^2 _ ^ 0/2 
4sinn7r2"n! 


<•(— i+,i— ) , 

x «■ - 1 m- zr-'-'dt. 

•'c 


(24) 


The functions T a n and V n differ slightly from those defined originally 
by Gegenbauer. They satisfy the equation 

{i - 1 ) w " + 2 (a -i l)zw' - n(n + 2a 4- l)w = 0. (25) 

But (25) is satisfied likewise by each of the following eight functions: 

ft , , (* 2 - lPftW , (S 2 - ir^n-l I 

x a n , , (* 2 - mw* , (* 2 - 1 reu • 

Since the second order equation (25) admits only two independent solu- 
tions, certain linear relations must exist between the eight functions 
above. These follow directly from the theory of the hypergeometric 
function, and they are readily deduced from corresponding expressions 
for the Legendre functions given by Hobson. 

ft = 

ft*. - Ift- 1 

SC = e ^.. (” + 2a l ! (/ - i)-scn*. 

nl 


(26) 

(27) 

(28) 


jLLn— la — 1 ® 


Zati 


(-n - 1)1 
(—n — 2a — 1)! 


(* 2 - l)“ a 5tr*_i 


(29) 


y a = fL — [sin (n + 2a) ir ft - sin tot ftn- 2 a-il ( 3 °) 

n ir cos (to + a)* 

s E. W. Hobson, “Spheroidal and Ellipsoidal Harmonics,” pp. 188 and 195, 
Cambridge University Press, 1931. 
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m — a 

■t n+2a — 


— L BUi /wr 

it cos (n + o)7r 

ri (n + 2a) ! sin (n + 2c 
7rn! cos ( n + a)ir 


f 2 n -a T -a _ e sin nir fr _ 
^ ' n+2 ° w(n + 2a)! cos (n + u)tt 

2 — ari • /T»a “t - 2a) ! (2 j \—a rp 

-e sin airTn = T n i \ z 1 


i (n -j- 2a)x 3Z°„-i] ; 

(31) 

-KW - 2=°n-i] 

(32) 

[£n — 2^- n -2a~l] 

(33) 

irra. 

(34) 

1)! i), ( z ! - lriCU 

(35) 


- * - (_ „ " 2o - 'i ) i (2 - ^ (35) 

The factorial has been employed in place of the gamma function for 
non-integral as well as integral values of the argument, so that 
o! = T(a -f 1). 

In addition to these formulas connecting the various solutions of (25) 
frequent use will be made of the recurrence relations 

2 Ti = ■> . o K» + 1 ) r ”+ 1 + <« + 2 “> r ”- J (36) 

2w + 2a + 1 

{i - l) 4- n = (« + d n + 1 - (» + 2 a + i)z r» ( 37 ) 

az 

which apply also to . 

From (21) and (22) and the known properties of Legendre functions 
the solutions of (25) can be expressed in terms of hypergeometric func- 
tions. Thus for | 1 — z | < 2, 


r » = (=^- irf V» 

while for | 2 | > 1 

_ e axi 2 a \Ar ( n “b 2a) ! 

(n + a + l)!(2*)*-w*« 


7i — a, n + a+ ljl - a; 2~ )’ 


n 2a 2 n 2a 


1 , ,3 1\ 

-> n + a+ 2'I>) 


where 


(0 ~ D» ^ (a + s-1)! (fe + 8-1)! ,. (40) 

F(a,b,c,s) ( a _ 1)1 (6- l)! s! (c + s — 1)1 
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Through the transformations of the hypergeometric function and the 
relations (26)-(35) one may establish series expansions for the functions 
T and X for the entire complex domain of 2 . 

1 .4 Gegenbauer Functions for Integer Values of a. Spheroidal wave 
functions associated with a complete spheroidal surface must be periodic 
in the equatorial angle In this case a is a positive integer or zero 
and will be designated by the letter to. The relations between the 
functions defined in Sec. 1.3 then reduce as follows. For all values of n 

we have 


fl 4-2771 


TZ = 


tv*-* 

(41) 


, (42) 

{n + 2m)! ^ 
n ! 

(43) 

(n + 2m) ! ^ 
n ! 

(44) 

(n + 2m)! ^ , 

n ! 

(45) 

tan nx ^ 

(46) 


X 


Numerical values of these functions in various regions of the z-plane 
may be obtained from the following hypergeometric senes: 

For | z | > 1, 


_ ±m (-2) Tm v / rr(n±2m)! 




£n (n ± m + £) ! (2 z)" +1 

X ( 2 2 - 0 Tm p( 2 


2 + 2 71 + 1 • 3 


- i ;»*» + !; J) 


(-2) ±m Vk (w ± 2ro)! 

(n±TO + i)!(28)-*"^ 

'n rfc 2m + 2 n i 2 m + 1 


(47) 


X F 


(’ 


; n ± m + 


■i) 


For | 1 - z 1 < 2, 

_ (n±_2m)l F /'_ n * + 2ro + 1; 1 + w; , 

2 m ?i! to! V 1 ' 


(48) 
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sc: = h n[in + 2m) ] " 8 _ 27 ~ \} z 

"P 1 ( - 1 y (— n— m 4-r-l)Un + m + r)\(m-r-l)\ ( l_^jY 

xl- h V 2 / 

r—0 

A (-» + r- 1)1 (« + 2 m + r)l ^ /l - »Y 

^ r!(ro-fr)! \ 2 / 

+ ( _ ir (jL±M ! ( z + l)-" 

' n\ 

- (_n-m + r-l)!(« + m + rt! ff(m + r )(LzfYl 
x 2^ r!(r-fm)! \ 2 / J 


(49) 


wherein n ^ 0 and 


<r(aO = t + \ + 


(50) 


as 


1 2 

For j z | < 1 and Itn(z ) 7 ^ 0, 


x 


-> 1 — rrt 


©' 


„ / n + 2m 4- 1 n . 1 . ? A 
F \ 2 ~~ 2’ 2 ,Z ) 


(51) 


+ 


zF 


rz = 2 m cos 


l ^ v r r (5+??)!2f 

n + 2m - l\, 

2 A 


/n + 2?n -J- 2 1 — n # 3 # «\ 

V 2 )' 


F 


/n + 2m + l — n_ 1. z\ 

V 2 ’ 2 ’ 2’ ) 


4- 2 m+1 sin ~ . 


71 4- 2 m\. 

" 2 /' „ ( n 4- 2m 4- 2 1 -jug g?) 

^©?> ( 2 • 2 2 j 


(52) 


( 


Usually we shall have to deal with functions for whxch n as well as m 
is a positive integer. In this case some of the above senes break off 
and hence are finite for all finite values of z. Thus when both n and 


8 



ELLIPTIC AND SPHEROIDAL WAVE FUNCTIONS 


are positive integers the various solutions fall into two classes : a set of 
integral functions, and a set characterized by singularities at z — il 
as well as at infinity. 


TABLE I 


r: 


I -n-2m- 1 


4— n— 2m— 1 




c* 2 - 

(s 2 - 

o* 2 - 

(/- 


7i and m 


i r 
i r 
ir 

i) 


i m— m 'i 

n+2m 


positive integers 

no singularities in the z-plane other than 

at z = oo 

singularities at z = =t 1 and z — 


If n is a positive or negative integer XT becomes infinite when n ^ 
-2m - 1 in virtue of the factor (ti + 2m)! in the numerator. It is 
finite if n > -2m - 1. Likewise TZ is zero when -2m - 1 < n < 0. 
It is finite for all other values of n. 


1 .5 Gegenbauer Functions for a = - b Comparable in importance 
to the case a = m is that in which a = ~b We have seen that this 
leads to the Mathieu functions. If we let z = cosh \fs, Eq. (20) reduces 
for a = — | to 

Hi - 617 = 0 (53) 

dv 

satisfied by hyperbolic functions. From the definitions of TZ and X n 
it can be shown that 

^ = (54) 

(jf - i)*iU = i < 55 > 

= ]/l * (56) 

T~* = - A ft sinh (57) 

n Y TT 

(z - 1) 4 T;_1 = ^ft^ cosh n\J/ (58) 

for all values of n. The relations between the functions reduce to 
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r? = Tzi 

(59) 

S? = --(** - D ! £U 

n 

(60) 

r? = - - ( + azt) 

7T 

(61) 

{i - l )‘ tU = - (az 1 - 3=1) 

7 r 

. Pi "" 

(62) 

= £ z? - - - i )‘rU 

2 L n 

(63) 


These are all integral functions, having no singularities other than at 
2 = co. 

II. General Spheroidal Functions 

2.1 Definition. In Sec. 1.3 it was noted that Eq. (7) reduces to (19) 
when 2 » 1 and c — > 0, and is then satisfied by 

U = (cz)~‘~ i Zi +M .i(cz) (64) 

where Z p (cz) is any Bessel function of order p and argument cz . Like- 
wise when c — » 0 Eq. (7) reduces to (20) and is satisfied by T 1 or 
provided b has the characteristic value 1(1 -f- 2a -f- 1). In general we 
shall express b in the form 

6 = Z(Z + 2a + 1) + €j(a, c) (65) 

where ei is a function of the parameters l, a and c, and vanishes with c. 
A circle draw r n about the origin as a center and passing through the 
singular points at z = ±1 divides the 2 -plane into two domains and 
limits the convergence of most series representations. We consider two 
systems of solutions. Those solutions of (7) that are valid in the neigh- 
borhood of the essential singularity at infinity we designate by U. They 
shall be represented in terms of the Bessel functions (64). Those solu- 
tions about the ordinary point z = 0 and converging at least within 
a circle of unit radius are to be designated by V, and will be repre- 
sented by expansions in Gegenbauer functions. These representations 
are remarkable in that only one set of expansion coefficients need be 
determined. The problems of numerical computations and of deter- 
mining the connections linking the solutions in one domain with those 
in another are thus vastly simplified. 
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2.2 U-Functions. The solutions of (7) valid when z » 1 will be 
defined by the expansion 

U(a, c; z) = (cz) _a_} £ a„Z n+p+a+t (cz), (66) 

where n is an integer and 0 <! p < 1. Upon introducing this series 
into (7) and applying the recurrence relations for Bessel functions it is 
found that the coefficients a n satisfy the three term recursion formula 

(w- + P 4~ l)(n + P + 2) 

(2 n + 2p + 2a + 3) (2 n + 2p -f 2a + 5) +2 

. (n + p + 2a - l)(n - j~ p 4~ 2a) 

(2rc + 2p + 2a - 3)(2ra + 2p + 2a - 1) 

b — (n p)(n p + 2a -j- 1) 
c 2 

2(n + p) + 2(ti -f- p)(2a -f- 1) -f- 2 a — l"l 

” (2n + 2p -b 2a - l)(2w + 2p + 2a + 3) J a " 

When c — » 0 there remains in the limit only 

[b — (n + p)(n + p + 2a + l)]a„ = 0. 

If then we choose 

b = (l - f- p)(Z + p + 2a + 1), 
where l is either an even or an odd integer, it is apparent that 

a n — > 0, for n H l 

and U reduces to a single term of order l + P + a + When c differs 
from zero the summation extends from — °o to + °o . As the subscripts 
of o n in (67) differ always by an even number, the even and odd series 
are distinct: the summation extends over even values of n when l is 
even, over odd values when l is odd. 

In order that (66) shall be an analytic representation of a solution it 
is necessary that the series converge for cz > 1. The maximum absolute 
magnitude of a Bessel function is independent of its order when cz — > = 0 . 
The coefficients a n must therefore converge in both directions from ai , 
the dominant term. This condition, in turn, restricts the separation 
constant b to a set of characteristic values bi . 

The quantity p plays a role analogous to the roots of an indicial equa- 
tion in the power series solution of a differential equation. By proper 



(67) 

= 0. 
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choice of p it is possible to terminate a certain number of the expansions 
at 7 i = 0 and thus to restrict the summation to positive values of n. 
However, this is untrue in other cases and in general the summation 
must be extended from - » to + «> , over either even or odd values of n. 
An earlier paper 6, failed to take this into account and the conclusions 
drawn there must be modified accordingly. 

If now we let 

U(a, e; z) = (s 2 - 1 r a X(a, c, z), (68) 

it is apparent that X must satisfy the equation 

(z 2 - l)X" + 2(1 - a)zX' + (cV - b - 2 a)X = 0 (69) 

and X can be represented by expansions of the form 

X{a, c; z) = (<*)“"* £ (&- Z n ^ +i {cz), (70) 

n ' 


where the a' n > satisfy the recursion formula 

(n' + l)(n' -f 2) 0 ' 

(2n' - 2a + 3) (2 n' - 2a + 5) n + 

(n' — 2a — I) (w' — 2a) 


+ 


+ 


(2 ft' - 2a - 3) (2ft' - 2a - 1) 
b + 2a — (ft' — 2a + l)ft' 


a n '_2 


(71) 


i 


c 2 


2 ft' 2 + 2ft' (1 - 2a) - 2a - 


T)] al ' = a 


(2ft' - 2a -f 3)(2ft' - 2a - 1). 

Upon comparing (71) with the recursion formula (67) we observe that 
the two are identical if 


ft' = ft -f" p -f* 2a, 

(ft + p) ’• 


ci n i — 


(n 4- p 4- 2a)! 


Ca n , 


(72) 


where C is any constant. Since the Bessel functions of negative order 
satisfy the same relations as those of positive order, further solutions 
of (7) are obtained involving the same coefficients a n . 

We shall define a set of eight spheroidal functions as follows in terms 
of the Bessel functions of the first and second kinds. The normaliza- 
tion of the coefficients will be discussed later. 

Ui{a, c;z) = (cz) _ ° _i ]T/ anJ n+ P +a+i(cz) (73) 


6 J, A. Stratton, Proc. Nat. Acad. Sci., 6 , 51—56, 1935 and 6 , 316-321, 1935. 
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Ulifl’t C‘, Z) — (cz) ° * E ClnJ—n—p-a~\{ CZ ) 

» 

U 3 (a, c; z) = (cz)“ *(z 2 —1) E p _j_ 2a) ! 

^ c ?-») = (c2) “" i(z2 “ 1)_a ?' (^+7T^! 

U 6 (a, c; z) = (cz)'*"* E' a«iVn + p+a+i(cz) 


17, (a, c; z) = (cz) a * E' a n N -^p- a -\(cz) 

n 


l/ 7 (a, c; z) = (cz) a *(z 2 - 
*7 8 (a, c; z) = (cz) a_i (z 2 - 


l)" a E' a. 

n 

ir a E' ». 


( n 4 p) ! - 
(n + P + 2a) ! 

(n 4 p) • 

(n + p 4 2a) ! 



(74) 

J n+p+a+\(fiZ) 

(75) 


(76) 


(77) 


(78) 

N n -f-p4-a+i(cz) 

(79) 

N-n-p-o-iW- 

(80) 


The coefficients a„ of all these solutions satisfy the recursion formula 
(67), and summation is extended from negative infinity to positive 
infinity. The prime over the summation sign indicates that only even 
values of n are to be taken if the index l of the separation constant is 
even, only odd values if l is odd. 


2.3 T -Functions. We construct next a set of solutions for the do- 
main | 3 | < 1. When c = 0 Eq. (7) is satisfied by T l+P {z) provided that 
the characteristic value b t belongs to the set (l + p)(l + P + 2a 4 1). 
If, therefore, c does not vanish, one is led to try solutions of the form 

F(a,c;z) = dnTUpW, C* 1 ) 

n 

where n is an odd or even integer and 0 ^ p < 1. The I th term of the 
series is dominant when c is small. Upon substitution of V into (7) 
and applying the differential and recurrence relations of the Gegenbauer 
functions it is found that the coefficients d n satisfy the recursion formula 


(n -f P 4 2a 4 2)(n 4 P 4 2a 4 D d 
(2n 4 2p 4 2a 4 5)(2 n 4 2p 4 2a 4 3) n+ 


4 


4 


4 


d n . 


(n 4 p)(n + p — 1) 

(2 n + 2p + 2a - l)(2n 4 2p 4 2a - 3) 

f 2(n + p) 2 + 2 (n + p)(2a + 1) + 2a - 1 
1_ (2 n 4 2p 4 2a — l)(2n 4 2p 4 2a 4 3) 

(n 4- p) (n 4 p 4- 2a 4 1) — b ~| ^ _ q 

& ~ J 


(82) 
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By comparison of (82) and (67) it is easily shown that, apart from a 
proportionality constant to be determined by the normalization, 


d n = i n 


(n + p) ! 

(n + p + 2 a ) ! 


(83) 


The recursion formula (82) is unmodified when n + p is replaced by 
— n — p — 2a— 1, so that two more solutions can be constructed from 
the functions T“ n _p_ 2a _i , X- n - p ~ 2 a-i • Furthermore, as in the case of 
the [/-functions, another four solutions are obtained by the trans- 
formation 


F(a,c;z) = (z 2 - 1 (84) 


where Y (a, c; z) satisfies (69). Thus eight formal solutions of (7) are 
established which we designate as follows: 


VM c; z) = (** - l)-“ £' (85) 

n 

y 2 (a, c; *) = ( s 2 - l)-“ E' a n r Strife) (86) 

n 

7, (a, c; z) = E' aS ^ + , p) J (z) (87) 

7Aa,' i .)-Z'a.C J £±^XL tt W ( 88 ) 

F 6 (a, c; I !) = (*•- 1)- E' (89) 

n 

7, (a, e; *) = (z 2 - 1)““ E' *»*" T^ +p+ ^(z) (90) 

n 

v ’<* * •> - ?' ^ W (91) 

(92) 

Again the prime over the summation sign indicates that odd or even 
integer values of n are to be taken according as l is odd or even. Unless 
otherwise specified the summation is from negative to positive infinity. 

For large values of | n | the absolute magnitude of the ratio T a n /T a n t . 2 
or X a n /X a n +2 approaches unity at either z = 0 or z = =bl. Thus the 
convergence of the F-functions, as was that of the [/-functions, is deter- 

QO 

mined by the convergence of the sum a„ . 
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2.4 Calculation of the Coefficients a n and the Separation Constant: 

For sufficiently small values of c one may express the separation con- 
stant b and the coefficients a n as power series in c. These series are 
then substituted in the recursion formula (67). Upon collecting terms 
and equating coefficients of like powers of c to zero one obtains in a 
single process both the characteristic values b and the expansion coeffi- 
cients a n . However, the convergence of these expressions is usually 
slow and they prove unsatisfactory for values of c 2 > 1. At present the 
most useful treatment of three term recursion formulas is based on a 
method of continued fractions. 7 W e shall describe briefly the procedure. 

To avoid cumbersome factors let us write the fundamental recursion 
formula (67) in the form 


4n+ 2tt n +2 T B n d n + C n -2«n-2 — 0 


where 


. (n + p ) (n 4- p — 1) 

An ~ (2 n + 2p + 2a - l)(2n + 2p 4 2a + f) ’ 


B n = 


b — (n + p)(w + p + 2a -fl) 


2(n + p) 2 4 2(n 4 p)(2a 4 1) + 2a 1 
(2 n 4 2p 4 2a — 1) (2 n 4 2p 4 2a 4 3) 

„ (n 4 p 4 2a 4 l)( n 4 p 4 2a 4 2) 

Cn ~ (2n 4 2p 4 2a -f l)(?n 4 2*> 4 2a 4 3) ' 

Evidently (93) can be developed as a continued fraction. Thus the 
ratio ai/ai-i is 


B, + A, +2 ^t? 
Ci l — 2 

A1+2C1 


Bi + 2 — v 


-A 7i Cn— 2 

jB n 4 A n+ 2 


v E. L. Ince, Phil. Mag., 547-558, 1928. 
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If the series is to converge, the value of ai/ai-i will be unaffected by 
the last term a n +i/a n provided n is chosen sufficiently large. In practice 
one may usually start with n ten or twelve integers larger than l, neglect 
<»{+i 2 /aj+io > and upon running up the fraction an expression for the 
initial ratio ai/ai-i is obtained in terms of b which is contained explicitly 
in B. Repetition of this process with a neighboring value of n will 
indicate whether or not a sufficient number of terms have been taken 
into account. 

Next the ratio atjai-i is approached in the opposite direction. Eq. 
(93) can be developed in the form 


«t-2 

di 


Ai 


0 , 1-1 


At 


Bi-i — 


A 1-1 C i -4 
Bi-i — \ 


A— n A.2C— n 


B- n + C- n -1 


fl-n-2 

a-n 


(96) 


Again the ratio a_ n _ 2 /<i_„ can be neglected for a sufficiently large value 
of n, and upon running up the fraction an expression is obtained for 
ai-i/ai in terms of b. The reciprocal of this ratio must equal the value 
of ai/ai-i obtained from (95). Thus (95) and (96), when properly re- 
lated, constitute a transcendental equation in b whose roots are the 
characteristic values. When these values have been determined, the 
continued fractions provide an easy means of computing the coefficients 
a n in terms of any arbitrary coefficient a t . 


2.5 Termination of the Series. In general the series expansions of 
the U- and F-functions extend from minus to plus infinity. However, 
by a proper choice of p they may be terminated in certain cases at 
n = 0 or n = 1 and extended over positive values of n only. Thus for 
the even series 


CL— 2 
Clo 


Ao 

B - 2 -f- C 


CL — 4 
0,-1 


(97) 
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with 

i - ZiA - (98) 

0 (2 p + 2 a - 1)(2 p + 2a + 1) 

Then if ao is finite and a j* a_ 2 vanishes as p -> 0. Moreover, it 
follows from (96) that all the remaining coefficients a„ with negative n 
vanish in the same manner, while the ratio a„_ 2 / a n remains finite. 
Since o_ 2 = 0, the coefficients with positive subscripts can all be ex- 
pressed in terms of Oc . From (96) and (95) one obtains 


02 _ Bo 
cio At 


a* 

Oo 


Co 


b, + aA' 

at 


(99) 


As in Sec. 2.4 the separation constant b is to be determined such that 
the two expressions for the ratio at/ ao are equal. The treatment of the 
odd series is identical. 

Since the functions J v {cz ) and N p (cz) are finite for negative orders p, 
it is apparent that the functions Ui, U 2 , U 6 , and U 6 defined in Sec. 2.2 
are summed only over positive values of n, starting at n = 0 or n = 1 
as l is even or odd. The remaining functions U 3 , Ui, U 7 , U 8 con- 
tain terms involving negative values of n in virtue of the factor 
( n + p )!/(n 4- p 4- 2a)!. Thus if a = m, a positive integer, the limit 
of a n (n + p)l/(n + p-f 2 m)!asp-* 0 is finite when n > -2m - 1, 
but zero for all values equal to or less than this number. The series 
again break off on the negative side, but at n = —2m and n = —2m 4- 1 
for the even and odd cases respectively. When a is a negative integer 
the series also terminate. If a = -m, the even and odd series start 
at n — 2m and n = 2m + 1 respectively, where m is a positive integer, 
and extend to infinity in the positive direction. Finally, when a is not 
an integer, all terms of the summation, both positive and negative, must 
be included. 

Similar considerations hold for the F-f unctions. If a = m, the func- 
tions Vi, V*, V bt V ti V 7 , V s are summed over positive values of n 
only, starting at n — 0 or n = 1. Actually Vx and V 3 are infinite in 
virtue of a factor 1 /p, but since this factor is common to each term of 
the series it can be removed. The summations for the functions V 2 
and Vi extend from negative to positive infinity, even when a is a positive 
integer. The ratio a n /p for negative values of n was shown to be finite. 


17 


L. J. CHU AND J. A. STRATTON 

Since T% P+ U*) for n < 0 and ££+, for n < -2m contain a factor 1/p, 
the products 

m ! "f P) • rpm 

a„SCr P+! ». and ( V+2 m + p) ! " + ' 

are likewise finite. 

2.6 Mathieu Functions. The special case a = -£is of particular 
interest since it pertains to the Mathieu functions. Taking p = 0, one 
now obtains 


A n — C n — 



and the recursion formula reduces to 

flM1 + o^-[a"+J(rf-6)]».-o. (ioo) 

Unlike the other cases, this relation is the same for both positive and 
negative values of n. Moreover, the Bessel and Gegenbauer functions 
are now either even or odd functions of the independent variable z. We 
shall construct V- and F-funetions that satisfy one or the other of two 
conditions. For the first type we shall take 

' a„ = (— l) B a_». (101) 

The series fails to break off, but in virtue of the symmetry of positive 
and negative terms it can be summed from n = 0 or n - 1 to n - 
In this case Oo differs from zero. When the index l of the separation 
constant is even one obtains b from the relation 

a* i 26 . 1 _ (102) 

do c 2 ABi ■+■ n */ a * 

which can be expanded as a continuous fraction as in (95). Likewise 
when l is odd, the characteristic values of b are given by the roots of 

- 3 = 3 + 

Ol 

For the second type we 


?(* - » = 


4B3 -b ns/ Os 

shall have 

On = — (— l)”o_n • , 


(103) 


(104) 
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Obviously a 0 = 0 and the characteristic values of b corresponding to 
even values of l are obtained from 


°- 4 = 2 + i ( 4 - 6 ) 
a 2 c 2 

while the roots of 


1 

4B 4 + cl(,/cl 4 


(105) 


a* 

ai 


1+-,(1 - b) = 
c 2 


1 

4B$ + ctt/aa 


( 106 ) 


are the characteristic values of b associated with odd values of l. The 
relation of the functions constructed in this way to the Mathieu func- 
tions defined by other authors will be shown below. 

The U- and 7-functions themselves have certain properties that are 
peculiar to the special case a — — f, p — > 0. We now find 
(n -f p)!/(n + 2a + p)I = n, n + p + a + \ - n, and since for any 
Bessel function of integral order Z n = ( — 1 ) n Z- n , there are two cases 
to be considered. 

1. Even Functions a n — ( — 1 ) n a- n .‘ Then U 8 , U 4 , U 7 , U 8 vanish, 
for the positive half of the series is equal but opposite in sign to the 
negative. The remaining functions U\ , C/2, C/$ , Us are of the form 

X)' a n Z ±n (cz). 


Likewise V 7 and 7 8 vanish, while the remaining 7-functions can be 
put into the form 

X/ a n i n cosh nt , 

n 

where t = cosh -1 z, and are therefore proportional to each other. Thus 
when p = Owe have only one independent 7-function. A second solu- 
tion in this case must be found by a limiting process.. Mathieu func- 
tions constructed in terms of cosh n\p alone are even functions of and 
it is apparent that this has nothing to do with the odd or even character 
of l. 


2. Odd Functions , a„ = — ( — l) n a_ n .• The functions U\, U 2 , Us, 
U 8 now vanish while C/3 , C/4 / U 7 , C/ 8 are all of the form 


V z 2 - 1 
cz 


^ ^ TlClnZ^n’ 

n 
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Likewise F 6 and F 6 vanish and the remaining F-f unctions are all of 
the form 

52' i sink n $ 

« 

and hence proportional to each other. They are odd functions of \f/, 
for even as well as odd values of l. 


2.7 Convergence. The recursion formulas associated with equa- 
tions of Mathieu type are of three terms. Consequently the conver- 
gence theory is more difficult than that of hypergeometric equations 
whose recursion formulas are always of two terms. In the present 
section we shall indicate the nature of the convergence and the conditions 
which appear to justify certain operations carried out in the remainder 
of this paper. The correctness of the results is supported also by ex- 
tensive numerical calculations. A more thorough investigation is 
nevertheless desirable. 

According to (93) the recursion formula relating the expansion co- 
efficients can be written in the form 


a n 


C 


Gn-2 


Bn + A„+, 


(107) 


Let r n = a„/a n -2 . Then as n — ► 

c 2 

r " 2c 2 - 46 + 4n 2 - cV n+2 * 


(108) 


There are two possibilities. The separation constant b may be such 
that the denominator of (107) approaches zero with increasing n. The 
ratio r n increases as n and the series $2 a * diverges for all values of c. 
But b may also be chosen such that r„ +2 0. This was the procedure 

described in the foregoing sections and it is now apparent that for these 
characteristic values of 6 



(109) 


The series of coefficients a n converges accordingly. 

It is interesting to note the origin of these two sets of b, the one 
leading to divergent, the other to convergent series. The recursion 
formulas for the coefficients constitute a system of difference equations. 
A two-term difference equation is the counterpart of a first order dif- 
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ferential equation. One is free to specify the behavior of the solution 
in the neighborhood of n = 0, but the character of the solution as n — » °o 
is then fixed. A three-term difference equation corresponds on the other 
hand to a differential equation of second order, to which there are two 
independent solutions. Thus we may specify the value of the ratio r n 
for small values of n, but this condition is no longer sufficient to deter- 
mine the behavior of r n at infinity. From the theory of difference equa- 
tions it is easy to find asymptotic solutions when n — > <» . The reader 
may verify that there are in fact two, the one becoming infinite, the other 
vanishing with increasing n. Of these two we choose in effect the latter, 
continue it analytically into the region of small n, and then select b 
such that the function r n assumes the proper value at n = 0 or I.. 

Consider now the U- and 7-functions and examine the behavior of 
Bessel and Gegenbauer functions whose order is very much larger than 
the argument. Let cz — p sech a 5= p. As p > 00 •, one must have 
a — > °o if cz is to remain finite, or tanh a — » I. Then according to 
Watson 8 , when a is any fixed positive number and p is large and posi- 
tive, one obtains the following asymptotic representation: 

p(tanh a— a) 

J p (p sech a) — » = ==■ . (HO) 

•\/ 27 rp tanh a 

Since a » tanh a, (110) vanishes exponentially with increasing p. On 
the other hand 

p(a— tanh a) 

N p (p sech a) — > 6 , == • (HI) 

V hirp tanh a 


This function grows large without limit as p —* and the same is true 
of J- P (z ) when p is not an integer. Nonetheless the corresponding 
[/-functions converge in virtue of the rapid decrease in the coefficients 
o n . Thus ifp = n + p + a+ 5 , 

a n Np(cz) ^ ^2 (a— tanh a) ( 112 ) 

a n -zN p -i(cz) 4n 2 


as n — » oo. The asymptotic formulas (110) and (111) hold only when 
p » cz » 1. Hence the convergence of the [/-functions has been 
demonstrated only for large values of the argument. Presumably the 
expansions (73) -(80) are convergent throughout the entire domain ex- 

8 G. N. Watson, "Theory of Bessel Functions,” p. 243, Cambridge University 
Press, 1922. 
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, 'mn z = 0 From a practical 
tenor to a circle of unit radius about the onpn ^ ^ ^ only when 

point of view 1 however t fl*» whe n the order p is equal 

z is large. The functions N,(cz) gr firgt few terms increase 

’J & - 01 " 

This is easily shown when cz <5C 

(113) 


a n N P (cz)_ _ 
a n -iN h ( w ) 


1 

— - 


, __ w nr a = lihas beeI1 sb ° Wn ^ 

Exceptions occur when a - ted. All expansions in terms 

in these cases certain . senes^an ^™ted ^ ^ only 

of Bessel functions of th including the points z = ±1. 

converge over the entire *-ph ane ndud^ fay ^ behav i 0 r of 

The convergence of the T-fu ^ required asymptotic repre- 

Gegenbauer functions of large # <* ■ coefficients o, diminish accord- 

sentations are given 1 by , H b °u at the expansions (85-92) are convergent 
ing to (109), it can he shown tha t the < ^ ^ ^ origin and in 

both in the intenor of a circle £ ^ ^ z = ±1 the value of 

the region exterior to that ^ ^ to a logarithmic factor or a 

the functions may become minute owing 

factor of the form (z - 1) • 

2.8 Expansions of V- T^ytc' coined 

mains the fundamental n^nthe prtSectioL. This will be 
between the 16 functions e n various functions m terms of 

accomplished by first representing thev va ^ ^ ^ , , | > l 

power series; coefficients can ? { Besse i a nd Gegenbauer 

we shall employ power sen^ -P"t Terms may be grouped and 
functions that are absolute^ ccnverg ^ ^ the coefficients 

the order changed arbitrarily. shown that the power senes expan- 

o„ form a convergent series . this wa y likewise converge, 

sions of U- and "f that in general only asymptotic 

^ Hobson, loc. cit., PP- 302-313. 
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externally by a circle of very large radius excluding the point at 
infinity. 

Within this annular region we find for the first four U- functions: 


”-r-‘ (?)'?■' ■■(!)' 


xE 


^ k\(n + p + a + k + ?) 


-,(¥)■ 


(114) 


Ih = 2" a_i ' C - 


(f 


xE 


1 


ICZ 


—a <> a I 


u* = ( * 2 - ir a 2 


fc =o k\(—n — p — a-|-fc — |)!\2 
(n -f p) 1 / czV 


2A: 


(115) 


/«Y° + ' y- 

(2) ? a ” 


x Z 


(n + p + 2a) ! \ 2 
1 


lh = (z 2 - 1)"°2' 


—a o 0— i 


^ofc!(n + p + a + fc-t-^ 
(n 4- p) I / cz 


*)! V 2 / 


(116) 


dr?'- 


(7i + P + 2a)!\2 


xZ 


1 


rofc!( — n — p — a -b & — f)!\2 


jcz 


2k 


(117) 


When | ^ | > 1 the ^-functions are represented by the series (39). In 
this domain the first four 7-f unctions assume the form 

Vx = (-iP V~A2zY E' <*. 2 a - -it ! (a,)’ 


(— n ~ P - 1) ! 

x V' (2r - n - p - 1) [ Qg) 

~o r\{r - a - n - p - §) ! 


(118) 


— 2r 


— 2r 


* = ( -i)-vi(2,)— j E- a. I-*)" 

y (2r + » + p + 2a)! (2z) 

r =o r! (r -H n + p + a + 2 ) • 

Fs = (-*rvW - i)-( 2a )^ sin(2 ° + ,>) - E' «•(«*>’ 

SUl p7r n 

y (2r - n - p - 2a - 1) ! (2z) -* 
r “o r ! (r — n — p — a — I) ! 


(119) 


(120) 
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Vi = (-r Vi rfe 2 - l) _a (2z ) _P_1 E' 0.(-»2*r 

fl 

(2 r + n + p) ! 


xL 


TTi 


— 2r 


(121) 


. ^ r ! (r + n + p + a- + 2 ) • 

Let s be an integer which differs from n always by an even integer. 
The eight solutions can be arranged as follows: 


_ o— » ( 


U i = 2 


V/ / CZ 


f) ?\a 


x E' 


a n i 


( 122 ) 






s — n\,f s n l 


+ P + a ) • 


x r , , 

ng« ( n — s 


dj\% 


(123) 


s n + 1 


- P-a y. 


Ui = 2 -V - 1) 

xE' 


-fpE'd 


a„i* n {n + p) ! 


(124) 


(i^y (i±±±l + a + p)l (» + p + 2a) ! 


u..rV-o-(tTs(g)' 


xE' 


a„i"-(n. + /.)! 


(125) 


n >» ( n — ( s -f n + 




2 y\ 2 

V, = (-|)“ a \/7r (22 ) P L" ( 2z ) 8 ( _s - p - I)’- 


(n + p + 2a) ! 


xE' 


a n i n ( — n — p — 2a — 1) ! 


(126) 


s + ” +1 - a- p)u-» - P- 1)! 

y 2 = (_§)-“ V^ (22)“ !o “ 1 "' > E' (2*r(« + 2o + #>)! 


xE' 


a n i n (n T p)! 


(127) 


"*• ^ s + n+ 2 4. a + p)> (n + 2a + p)! 
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X E' (2z)*(~ S - 2a - p - D! 


( 128 ) 


x E' 


a n i 


ni . _ a - p ). 


y, = (-J)V^ (2* - ir*(2*r' _1 E' (2»)"(« + p)1 


XE'/ 


a„i 


(129) 


s-n\ t (8 + u±_l + q + p 


’) 


2.9 Relations between the Functions. Since the eight functions 
defined in Sec. 2.8 satisfy one and the same equation, there can be at 
most two that are linearly independent within a common region of 
convergence. One observes that, so far as powers of z are concerned, 
each {/-function is identical with the 7-function of the same subscript, 
and we conclude that the two are proportional. Thus 

Ui = KiVi , t/ 2 = X 2 7 2 , (130) 

u z = K 3 Vz , Ua = KaVa . 

The proportionality constant K is called the joining factor.. Upon 
equating the power series expansions of U and V the factor K is given 
as the ratio of coefficients of corresponding powers of z. 

To obtain further relations between the various solutions we make 
use of the formula defining the Bessel function of the second kind, 


N p (x) 


Jp(x ) COS PIT — J-p(x) (131) 

sin pi r 

together with the functional relations between the Gegenbauer functions 
given in Sec. 1.3. Through these it is easy to establish the following 

connections. 


Relations among the U -functions 
U 5 sin (a -f p)x + (— 1)” ^6 


V i = 


u> = 


cos (a T- p)ir 
Un sin (a -H p)ir + (~~l) n ^ 


cos (a + p)tt 


(132) 

(133) 
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Uz = 

u< = 

Ui = 
Uz = 

u, = 

U 8 = 


Ui sin (a + p)^ + ( — 1)" C/ s 
cos (a + p)ir 

Us sin (a + pK + ( — 1)" C/7 
cos (a + p)r 

C/ x sin (a + p)t + (— 1)" C/g 
cos (a -f- p)t 

C/ 2 sin (a -j- p)* + ( — 1)" C/ i 

COS (d + p) 7T 

C/ 3 sin (a + p)t + (~l) n C/ 4 
cos (a -f p)t 

C/< sin (a + pV + (— l) n C/ 3 


COS (a + p)7T 

Relations among the V -functions 


-2a 


sm px 


7 S 


^ ^ sin (2a + p)ir 

7 2 = (-l)- 2a 7 4 

7 6 = 7e 

7 7 = 7 8 


t i\a 2 . T r SinpTT 

(— 1) -smaTKi = - — 7= — j — c- 
x sm (2a + p)jt 


7 7 - 7, 


(134) 

(135) 

(136) 

(137) 

(138) 

(139) 


(140) 

(141) 

(142) 

(143) 

(144) 


(- l) a - sin airV 2 -7,-7, (145) 

7 r 

(— l) a 7r cos (a + p)tt 7 6 = sin pir(7 4 — 7 3 ) (146) 

(— l)°x cos (a -f p)tt 77 = sin (2a + p)tt(7 2 — 7i) (147) 

Relations between XJ~ and V-functions 

The power series expansions of the functions Ui , C/ 2 , C/3, C/ 4 , 7i , 
7 2 , 73 , 7i converge within the domain 1 < | z | < «> when b has a 
characteristic value ensuring the absolute convergence of the coefficients 
a n . Since U- and 7-functions of equal subscript have power series of 
identical form, they must be proportional to each other and the propor- 
tionality factor Kj = C/,/7,- , (j = 1, 2, 3, 4), must equal the ratio of 
the coefficients of 2 * in the corresponding series. Thus there are an 
infinite number of expressions for the joining factors, for s is any posi- 
tive or negative integer differing from n or l by an even integer. 
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Ki = —= 


(-D 


a+tl 2 a+p 
c 


V2ir 2 2c,+a+p) ( — s — p — 1)! 


E' 


X 


^( i ^n) I ( S _±n±i + „ + p ) I 


(148) 


E' 

n 


l ( — ti — 2a — 1 — p)!a* 


s C-i-K 


s -f- ft + 1 


a - p )! (-7i — 1 — p) ! 


K, = 


(-D 


o — */2 «2 (*+a+l +p) 


•v/2tt c'+ 2a+l+p (s + 2 a + p) ! 


E' 


t"Ctn 


X 


- (v>(- 


s + 71 + 1 I . 

a pi! 


(149) 


E' 

»=?• 


i n (n p) ! a r 


( S -i-”) ! ( 1 ± f ±J + « + ') &* + 2a + p)l 


K* = 


(- 1 ) 


—o+a / 2 *+ 2 a+p 


V27T 2 2(,+0+/,) ( — s - 2a - 1 - p) ! 

™ *-*(— n - 2a — 1 — p)!a„ 

4- /o _ *A / 

i>?» 

X 


( 5 __^)! ^_+n+_l + „ + p y { _ n _ j _ p) , (150) 


E' 

nS« 


l n CL n 




/ i \ — a— »/2 2 2 («+a+l+p) 

y^/ 2 n (yi + p) 1 fl» 

4 (W)r(- s + "+ 1 - a - p)l(» + 2a + p)l (161) 

X -T— . 

1 u n 

Every value of s must, of course, lead to one and the same value of K, . 
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2.10 Summary of Relations when a = m. When a = m, a positive 
integer, and p -» 0, the 16 U- and F-functions fall into two groups of 
independent solutions. The following table is so arranged that all 
functions in a group or column are identical. 

TABLE II 


Group 1 

Group 2 

- P V i 

f 2 

— pF 3 

F< 

Fs 

Uz/Kz 

F e 

UJKk 

f 7 

(-i)’ + ” +1 U,/Kt 

f 8 

—pUi/Ki 

-pUt/Kz 

(—\y +m+l pUt/Ki 

(-iy +a+ ' P u a /K, 

(_1 )"+"+ l U,/K, 


By proper choice of s the joining factors reduce to the following. 
For l and n even: 


V— m — $ 


K, = - P ? 


ao 


K, = 


(w -f |) ! 

2 -m+i 


[E 8 (o)r 


K z = p 


cirri — 4) ! (p — 2m ) ! a_ 2 m 
2 m "*(m-!)!( P - 2m)!a_ 2m 


F 8 (0) 
[Fg(O)]- 1 


(152) 

(153) 

(154) 


where 


K t = 2’” +i (>ra + j) ! y, (0) 
7rcao 


y 8 (0) = y 8 (a, c; *) |._0 = -V “» / 2 ( • 

2V "”-° g + »)l 


(155) 


(156) 


For l and n odd: 


Ki = — p 


2~ m ~*ca 1 

(m + 1) ! 


[Fj(O)]- 1 


(157) 
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m+i 


k 2 = 


vim 


(158) 


cV-f)Kp+ 1 -2 m)!ax_ 2 m 

2 n ^dm - f)! (p + 1 - 2m)! m-j, (159) 

K S = ~ P — 7T 

(160) 


_ rVH)! r!(0)i 

ire 2 fli 


where 


Va(°) = ^ ^(a.cjz) 


*-o 2" 


2> 


IT 


Vt*- 1 


(161) 


The coefficients a_ 2m and ai_ 2 m vanish as p — > 0 and m > 0, but the 

products (p - 2m) !a_ 2 *» and (p + 1 - 2m )- a ^”> are finite ’ 

For both even and odd values of n 


KiKi = KA = ~ 


2p 

•7TC 


2 11 Summary of Relations when a - l- When a 2 and 

p - 0 the joining factors are reduced to the followmg by PJP- 
of *. These formulas apply both to the even and the odd Mathieu 

functions. 

For l and n even : 


= - i t / l oo [ ao + 2 S' a "] 

*•- - *■ \f\ [* + 2 S' a -] 

= i I f/l “ a-,) [S na ’l 

K t = - i| z^/ 2 (“• - [s n<L ’] 


For l and n odd: 


Ki = — 7= (°-i ~ °0 f 2 na "l 

2\/2t J 


(163) 

(164) 

(165) 

(166) 

(167) 
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K 2 = — 4= (a_! - at) 
2V2t r 


K 3 = — — -=. (<X-i + <*i) 

2V2^ 


-K4 = (fl-l + &l) 

2V2tt 


[s-r 
K-r 
rs-r 


(168) 

(169) 

(170) 


j Even Functions , a n = ( — l) n a_ n . The functions Us , U 4 , U 7 , Us, 
V 7 , Vz vanish as p — * 0, as do also the joining factors K s and Ka , and 

Kx = (-1 ) n K 2 . (171) 

The remaining functions again fall into two independent groups. 

TABLE III 


Group 1 

Group 2 

Vx 

u„ 

v 2 

(-1 ru. 

Vs 


-v 4 


. 7T T/ 

l ~ Vs 


2 




O 

Ux/Kx 


UsIK, 



As before, all functions occurring in one group or column are identical 
and represent one independent solution. The Gegenbauer functions 
for a = — §• can be expressed in terms of exponential and hyperbolic 
functions; by writing z = cosh ^ and doubling up the series one obtains: 
l and n even, 


Vx 


-VsD 


+ 2 a n i cosh nyp 

n=2 

00 

Ux = a 0 J 0 (cz) 

n=2 

U B = OoNoicz) rf 2 J2' a n N n (fiz) 




(172) 

(173) 

(174) 
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l and n odd, 


Vi = iV 2ir dni n cosh nj/ 

n— 1 

CO 

Ul = 2 Yf UnJnicz) 

1 

oo 

t/s = 2 E' OnJV„(c«). 


( 175 ) 

( 176 ) 

(177) 


TV 30 1 


V, and Cd are Mathieu ,. ;£ - ^h. J* ** £7*. 

SSTb <* iml .nd (177) » — d » to 

SZirhood of , = X been discussed by W*. but he fails togive 
the analytic connection with a corresponding solution i valid wh ® 3 . 

Tn the next section we shall obtain another form of the second solution 
and the desired connections through application of a limiting process to 
the yanking functions V, and 7. • p, , U. , 

V 5 ,V 6 , and the pining factors K x and K* vanish as p 

K 3 = (-1 )"K4. (178) 

The remaining functions fall into two independent groups as in Table 
IV, the functions in any one column being all identica . 

TABLE IV 


Group 1 

Group 2 

Vi 

u 7 

-V% 

(-1)^ 

v 3 


V 4 


i\Vi 


2 


ilVs 


2 


Us/Kz 


Ui /K, 



10 E. L. Ince, Proc. Edinburgh Math. Soe., pp. 2-13, 1914-P15. 
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Upon doubling up the series so that the summation extends only over 
positive values of n one obtains 

Vz = iV a n i n sinh n\p, (179) 

n=l,2 

Ua = 2 VIE! £' a n nj *(cs), (180) 

CZ n-1,2 

V, = 2 £' a „nl V.(e*). (181) 

CZ n— 1,2 

When l is even these series start from n = 2, and from n = 1 when Z 
is odd. As in the case of the even functions a second solution valid 
when | 2 | ~ 1 is missing, but will now be found by a limiting operation 
applied to the vanishing functions U 6 and U 6 • 

2.12 Second Solutions of the Mathieu Equation. The object of 
this section is to find a second solution of Mathieu’s equation in the 
domain | « | ~ 1 and to join it analytically with representations of the 
same function valid when | « | » 1. In Sec. 2.11 it was shown that 
y 7 and Us for the even functions, and Vz and V 6 for the odd functions, 
vanish as p — » 0. The zero is of the first order in p and hence can be 
removed if one divides by p before passing to the limit. Since p is an 
arbitrary constant, the functions obtained in this manner are still solu- 
tions of the differential equation. 

Even Functions , a n = ( — l) n o-n. Formally 


y 7 = U 8 -= a/- a n i n sinh (n + pty, (182) 

Y 7T „ 

where the a n are also functions of p. They vanish as p — > 0, since 
a n i n sinh nrp = -a- n i~ n sinh (-np). We now divide U 7 and U 8 by 
p, differentiate numerator and denominator with respect to p and pass 
to the limit p = 0. 

lim — = lim — = A ft- S' Hal sinh n\l> + a n yp cosh nf], (183) 

p-*0 p p — o p Y TT n 

where 


t 

a n 


d 



(184) 


In virtue of the relation 


a 


/ 

n 


-(-!)"<*-«, 


(185) 
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which will be established in the following section, the function defined 
by (183) is finite and does not vanish. It differs from all other E-func- 
tions in that it is an odd function of i {/, and it must, therefore, be a 
second solution of the Mathieu equation. It is this solution that was 
found by Ince 10 in another manner. 

To find a relation between U- and 7-functions of the second kind we 
make use of (136) for the case a = — 

Ui cos p7r — (-I) 71 U 2 


Ub = 


(186) 


sin pir 


Upon replacing Ui and U» by their equivalents in terms of 7-functions 
and allowing p — » 0, one obtains 



U, = - lim ~ [KiVi - 

TT p-* 0 dp 

(-1 YKiVil 

(187) 

Let 







K' = d * 

, V' 

_ dV 


(188) 


dp 

p— 0 

dp 

P“0 


Since Ki = 

(-l) n K 2 , 7i = V 2 , we have 




V t = - {[Ki - (- 

7 T 

-lYK'zWt 

+ KdV i - V’,]\. 

(189) 

From (85), 

(86), and (55), when p N 0, 





7i = H 

/IV- 

m^n+p)* 

(190) 



* e ~t,n+p)<P ' 

(191) 


Thus 


7j - Vi = W^tt Y! ? (<*» sinh n+ + a n ^ cosh nf] 


.7 7 . 7g 

= ITT = IT 

P P 


(192) 


It appears that the solution U 6 valid when z is large is in fact a linear 
combination of solutions of the first and second kinds established for 
the region | z | ~ 1. 

We proceed now to calculate the derivatives of the joining factors 
K. In (148) and (149) put a = and obtain 



L. J. CHU AND J. A. STRATTON 
«+p a n i a ~ n 



'2 c 
UA- r 

•3T \4/ 


S /s — 1 1 / s + ^ p 


K 1 = 


■) 


(-s - 1 - p)! S' /„ _ 


a„i"(n + p) 


L — (193) 


nSi 


(VX-H-*-')' 



—a—p 


'2 c 

7T\4/ n£« /n — s 


CL n V 


Kx = 




g»i n (w 4~ p) 

<- 1+ '>'5. (■-),(!+-"+,), 


(194) 


2 y \ 2 

For even values of n let s = 0. After simplification we obtain 

k[ •■= -k; 


= ^ In 4 


a 1 0.(1 — 2 — 

. 1 I £3 n 

-*yi 


(to 4“ 2 (195) 

n— 2 

2o« E' [«» + (l - ! ) + s}]| 

co “ 2 

do 4" 2 Z On 

n-=2 


where 


* 1 

C = 0.577216, <r(x) = Z “• 


n^l W 


For odd values of n let s = 1 in Ki and s — 1 in K% . Then 




-on (C + i) + 2 £ 


/ 0>7I 


n = 3 


n 2 - 1 


l 


+ CH Z 


flU 4- Ur 

/ 


Z' a « n 

n= 1 


1 4- wc 


71—1 


n=l 




(196) 
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In general 

K[ = -(- 1 YK, ( 197 ) 

U s = -K[Vi + iKi—. ( 198 ) 

X p 

Odd Functions, a n = — (— l)"a_ n . Formally 

F 5 = 7 6 = a/~ XT cosh ( n + pit- ( 1 Q 9) 

f 7T n 

They vanish as p — > 0, since a n i n cosh n\p = —a^ n i~ n cosh (—n{/). 
Upon division by p and passing to the limit p — * 0 one obtains 

lim — = lim — = A /- X' cosh nf + a n \f/ sinh n\{/]. (200) 

p— *0 P p-»0 P Y X n 

This is an even function of for 


(-l)V.n, 


( 201 ) 


as will be shown below. 

To establish a relation between U- and 7-functions use is made of 
(138). 

U 3 cos px — ( — 1 ) n U± 


Ui = 


sin px 


( 202 ) 


Replacing U s and Ut by their equivalent 7-functions and letting 
p — > 0, we have 


V-, = - lim £ [K, 7 s - (- lyKtV,], (203) 

x p -*o dp 


or 

[7, = - {[Ks - (-D’XllV, + KAVi - 7l]}. (204) 

X 

For arbitrary values of p 

t, = i J\ E' <*«*■“ • irW *, (205) 

74 = -i 4j/ 1 E' (206) 
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Thus 


Vl — V{ = i\Z2ir i n [o! n cosh n\p + a n \p sinh n\p] 


. r V 6 . F 6 
= n r lim — = n r Jim — . 

p— *o p p— *o p 


( 207 ) 


As before we find that C/ 7 , valid when z is large, is a linear combination 
of two independent solutions which have been established as useful in 
the neighborhood of z = 1 . 

The formulas for K 3 and K 4 reduce to the following when a = — J. 



K 4 = 



S' 


a n i n - a (n + p) 
(n - s\, / _s + n 

\ 2 A 2 



2 \/ 2* ( S + p)!E 



(209) 


When n is even let s = 2 for K 3 and s = — 2 for K 4 . There follows 
after simplification: 



( 210 ) 
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When n is odd let s = 1 for K 3 and s = — 1 for Ka . We obtain: 

. ((LcV + jf'Al 

= = K, In | + -ij ^ - 1 

4 ^ £' «n 

w «*“ 1 


«. t [ 

»-i L 

On 4- On | 

(" ; '). 


[.?>: 

2 


In general, 

K',= -(,-lYK' t , (212) 

U; = -K' a V, + iK s U (213) 

7T P 

2.13 Recursion Formulas ior a* when a = — £» P — * 0* Let' e = 
6 - c 2 /2. The recursion formula (100) for the Mathieu equation is then 

a n+ 2 + a n - 2 + ^ (e — n 2 )a n = 0. (214) 

c 

Or, upon placing 

J* n — d n /d n — 2 > (215) 

it may be expanded as in Sec. 2.6 in a continued fraction 
— c 2 /4 


€ “ n2 + 4 r "+ 2 


-c 2 /4 


e-(n + 2) i - e _ (Tt ; / 4)2 _.. ; 

Replace n by n 4- p, differentiate (216) with respect to p, and let p — » 0. 

= lim ( 2n - T r n-f2^ 

p^o dp c2 V 4 / (217) 

= — — r\ {2n 4* r 2 +2 [2(n + 2) 4- r 2 n+i (2(n T 4) 4- •••)]}• 

& 
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From (215) we obtain 




(218) 


The values of r n can be determined from (216), whence in virtue of 
(217) one may consider r' n to be a known quantity. To obtain each 
individual a n , one must know one particular a' n for a specified n, or 
establish an additional relationship between the a« . There are four 
cases: Even functions, l even or odd; Odd functions, l even or odd. 

Even functions, l and n even. According to (101) a n = {-\) n a^ n , 
and hence a n = a_ n . Consider the coefficients to be plotted as func- 
tions of the variable n + p. In view of the symmetry about the axis 
n p = 0, it is apparent that ah = 0, a» = —a - n . f Since in this case 
do differs from zero we learn from (218) that n 2 = cor 2 • The remaining 
coefficients a' n can then be found from the general formula 



Even functions , l and n odd. Then a n — n . The coefficient^ a n 
is an odd function of n and its derivative is therefore even. a n = a_„ . 
In particular a[ •= a'-\ , a\ = — a~i = — 1. From (218) we obtain 
a[ = —air'i/2 and the remaining a' n are found from (219). 

Odd functions, l and n odd. According to (104) a n = — ( — 1 )'"«-» j 
and hence a n = a-n , a n = -aL n . From (218) one finds a[ = W 2 
and the remaining coefficients again follow from (219). , 

Odd functions, l and n even. Then a n = —a- n , a n = a_ n , ao = 0. 
Upon differentiating (214) and placing n = 2 one obtains 


/ 2 e r 

fl 2 “ -^ a °- 

From (217) we find 

r J= -ln(4 + i{[8 + —I} = 

r 2 c 2 Co 

and likewise from (218) and (220) 


( 220 ) 


( 221 ) 


r l 

r 2 




a 0 



( 222 ) 


Since u 2 differs from zero one may multiply (222) by Co to obtain 

Go = OzK, (223) 
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a finite quantity. The remaining coefficients can now be determined 
from (219). 

In view of the symmetric or anti-symmetric relations between the 
coefficients the second solutions of the Mathieu equation reduce to the 
following: 

Even functions, l and n even , 


— = A /? \ Oa\p + 2 sinh C0S h 71 1) 

p y ttL n=2 

], (224) 

and for l and n odd, 


— = 2 d/- Z)' "(<*« sinh n f + a n f cosh n$). 

P Y T n-l 

(225) 

Odd functions, l and n even, 


— = j/? faj 2 yy i n (a' n cosh rnp + sinh mf) 

p V TTL *“2 

, (226) 

and for l and n odd, 


— = 2 i/? Ha' n cosh n\p + a n ^ sinh nf). 

p y % n“i 

(227) 

The summations are extended over positive values of n. These solu- 
tions are identical with those found by Ince 10 . However, Ince’s treat- 
ment expresses the coefficients a n and a n as analytic functions of c and 
its application is limited to small values of c. 


HI. Prolate Spheroidal Functions 

3.1 Functions of the First Kind. In the preceding sections, sixteen 
solutions have been found for the basic equation (7) and relations 
established that permit the analytic continuation of any one solution 
over the entire z-plane. We shall now specialize these functions for 
certain physical problems. Consider first wave functions of the prolate 
spheroid. In this case a = m, a positive integer, and the upper sign is 
chosen in Eqs. (3), (5), and (6). The F-functions are then functions 
of angle, or position on the surface of a spheroid, so that the variable is 
confined to the range -1 < z < 1. Radial variation is given by the 
C7-functions, with the variable limited to the range z >1. Those solu- 
tions that are finite at the poles z — rtl are termed functions of the first 
kind. Functions of the second kind are infinite at these points. 
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For the standard 7-function of the first kind we choose V 8 , multiplied 
by a factor i~ l . We define: 


7i?(c, z) = i l V& = 52 d l n T«(z), 

n= 0,1 


(228) 


where by (92) 


d l n = 


•n-j 


n! 


(n + 2m) ! 


a«. 


(229) 


The summation starts from n = 0 when l is even, from n = 1 when 
l is odd, and extends over alternate positive values of n. 

The coefficients will be normalized such that each spheroidal function 
reduces to the corresponding spherical function as c 0. It proves 
more convenient to carry out the normalization at z = 0 than at the 
poles z = ± 1. Thus it is desired that when l is even, 

7iV(c,0)= £' * n(0) = Tm, (230) 

71 = 0,1 


whence 

7iV(0, z) = 2TM, 


and when l is odd, 


d 

dz 


VL\\c, 2 ) 

1 2=0 


d 

dz 



For even values of l one has 



(231) 


(232) 


(233) 


(234) 
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Hence 


- fc 1 ) 

E ' ) -4 

n=0 ( Yl . \ i 

\2 m j 




(235) 


when l is even, and 


2 2m 1 + m ! 


(¥> 


when l is odd. This fixes the absolute magnitudes of the coefficients a„ 
The corresponding solution of the wave equation (1) will be denoted 

sL\\c, v). 

S^(c,v) = (l-vT n VSl{c,q) 

(237) 

= •£' di PZUv). 

71 “-Oil 

For the standard [/-function of the first kind we choose Ui with an 
appropriate normalization factor. 


USHe 


, z) = *i 


= Klicz)^ E' a ^'n+m(cz), 

n=0,l 


where /c* is the normalization factor, and 


jn+m(cZ ) — 7^/2” JrH-m+i( C2 )- 


The corresponding solution of the wave equation is then defined as 

*V(c,0 = if- lr^VSKe,*). (240) 

When c£ — » °o , 


. , ?~ l ■ ( t 1 + ™' 
j n+m {cv -> — sin ^ * 


so that 


i) - ^ - (« - ^ ') .§; 1 ” _ ‘ (242) 
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(Strictly (242) does not follow directly from (241), since (241) is valid 
only when the argument is very much larger than the order A more 
careful examination of the asymptotic behavior of aj n . UcQ on the 
grounds of Sec. 2.7 leads to the same result.) 1 he normalization factor 
will now be chosen such that for very large values of c$, 


, i / l + m \ 

Rml (c, £) -> £ sin ^ irj. 


Thus ki is determined by 


Kl = c 


S' r. 1 « n ~ l 

a n i . 

n<=0.1 


The functions Vfffe *) and USHfi, z) are different representations 
of one and the same solution. From the results of Sec. 2.9 the following 
relations may be established: 

For l even, _ 

, T m / x Kty/va o FiV(c, z) (245) 


T (X) , S _ Kl"V « r mi 

Uml ^ Z) ~ ^(m^FSTl yiV(c, 0) ’ 


and for l odd, 


C/iV(c, z) 


Kl\/ 7 r CQ\ V mi (c, z ) 


2 ,n+2 (m + 1) 




From the orthogonal properties of the Gegenbauer functions it may 
be shown that 


£ vSfc *)<i - & = £ ssfo ^(c, ^ dr} 


0, I * V 


_ A/ , 

2 n^.i ^ (n + 2 m)!( 2 n + 2m + 1)’ 


l = l'. 


3.2 Functions of the Second Kind. These are constructed from 
functions of the second group tabulated in Sec. 2.10. They are char- 
acterized by logarithmic singularities at the points z - ±1. We 
define : _ 


UZ\c 


«) = ** t\f\ 


= Klicz)-™ Un^n+fn( cz )j 

n=0.1 






'.r ‘ 
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where n n+m {.cz) is a spherical Bessel function of the second kind, 


n n+m (cz) 


]/ 2cz 


N n 4-7?i-t- i • 


The corresponding solution of the wave equation is then 

RS(c, & = (i? - i T' 2 u%(c, ()■ 

As c£ 


00 


(249) 

(250) 

(251) 

(252) 


rSIc, cos ( c * _ nr * 

It is desirable to define V${c, z) in such a way that 
UH}{c,z) _ U mlic, z) 

FiVM) nV(c,2) 

and to obtain for it a representation valid in the neighborhood of z = 1. 
From Sec. 2.10 one finds 

«SK«, *> « «. (253) 

This is combined with (252), and leads, when 1 is even, to 


ra>(e,*) = 


(4m 2 — 1) [y (V (c>0 )f7 j; 


or 


cao(p — 2 m)!a_ 2 m 
(2) _ Ki V 7T a 0 FS(c, z) 

C ' ml (C ’ 2) 2- +1 (w + 4) ! ViV (c, 0) ' 

Likewise, when l is odd, 

(47 

Jai(p + 1 - 2m) ! ai_ : 


(254) 

(255) 


yt*> (e d f j vSi (=, i)T V. , (256) 

W, Z) l rZnJn 4- 1 — 27 n)\ai- 2 m\_dz Jz=0 


and 


t/S(c, ^) = 


ki-\At 


ca i 

2 w+2 (m + |)I fd 




(257) 


The Wronskian of the two independent solutions for large values of 
cz reduces to 


USi^USS - US±V% = e-a - 2 "- 2 


(258) 
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By analogy with the Bessel functions one may construct solutions 
corresponding to the Hankel functions. 

Vmi (c, 2 ) = u22 (c, 2 ) + iUS (c, 2 ) (259) 

*) = UiV(e, *) - 2 ) (260) 

ViV(e, z) - riV(«, *) + »riV(c, 2 ) (261) 

ViV(e, z) = riV(c, 2) - 2). (262) 

Obviously the ratio of any 2 ) to the 7.,(e, 2 ) function of the same 

kind and order is constant and equal to 

K\y/ TT <*0 


2 " +1 (m + |)!Fi7(c, 0) 
Kiy/ ir cd\ 


ira+2 


{m 4 


L' 


vL'Hc, 2 )] 

3(3), 


for l even, 


for l odd. 


(263) 


z=0 


The corresponding wave functions are RiVfe 5) and ^’(c, i), whose 
asymptotic behavior when 

R%(c, *) 


p(4) 
■atm l 


(e,i) 



(264) 



1 

(265) 

ct 



If by f) and S ml ( c, ,) written without superscript one under- 

stands any solution of the four kinds, then the prolate spheroidal wave 
functions of Eq. (3) are 

If = B.i(«, {) S m i(c, v) «*”"'• (266) 

IV. Oblate Spheroidal Functions. 

4 1 Definition. Let VJifi, .) and U ml (c, z) be standard V- and 
[/-functions as defined in III. The kind of function will normally be 
indicated by a superscript. Then V m i(-ic, z) satisfies 

(1 _ Z 2 )V" — 2(m 4 l)zV f 4 (pi 4 cY ) V = 0, (267) 

while iz) is a solution of 

(z 2 4 1)17" 4- 2 (m + l)zU' - (bi - cV)U = 0. (268, 

These in turn are the equations of oblate spheroidal wave functions 
Eqs. (3)-(6). 
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,) = (1 - nY' 2 V ml {-ic, ,), (269) 

jM-w, *0 = ns 2 + l) m/2 !/•«(-«, <0, (270) 

TF = «»,(-*, *0 i?) e ±iw ". (271) 

In the oblate as well as the prolate case the normalization of the co- 
efficients at is such that in the equatorial plane transverse to the axis 
of revolution, r\ = cos 9 = 0, the spheroidal S m function reduces to 
the value of the corresponding spherical function. 

V. Mathieu Functions 

5.1 Mathieu Equations. In elliptic coordinates both the angular 
and radial parts of the wave function (9) satisfy an equation of the 
form 

(z 2 - 1 )U" + zU' + (cY - b)U = 0. (272) 

In the case of the angle functions z is confined to the range \z \ ^ 1, 
while for the radial functions z ^ 1. Upon making a change of inde- 
pendent variable, 

z = cos <p = cosh \f/, ^ — up, (273) 

one obtains 

+ (b - c 2 cos 2 v)V = 0, (274) 

difr 

- (b-c 1 cosh 2 +)U = 0. (275) 

d\f / 2 

Independent solutions of these equations will be defined below in a form 
which the authors believe to be the simplest and the most useful for 
application to physical problems. Apart from a slight change in the 
notation for the radial functions, these are the definitions proposed 
earlier by Stratton and Morse 6,11 . 

5.2 Even Functions of the First Kind. In this case the expansion 
coefficients are subject to the condition a„ = ( — l) n a_ n . Let 


D ° = *4/!“°’ 


O. = i” + V27rO„, (n > 0). 


11 P. M. Morse, Proc. Nat. Acad. Sci., 6, 66-62, 1935. 
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Then for both even and odd values of l Eqs. (172) and (173) reduce to 
Vi = E' D n cosh n\j/, (277) 

n— 0,1 

Ur = T r' , - l D„J n (cz). (278) 

V 7T n—0,1 

We take V x to be the standard V- or angle function of the first kind 
and define : 

&c f a) (c, cos <p) — 'jt,' D l n cos n<f>. (279) 

n— 0,1 

This function will be normalized at the point tp = 0, or z — 1, such that 

Se?\c, 1) = 1, (280) 

or 

E' Dn = 1. (281) 

n— 0.1 

For the standard U- or radial function of the first kind we choose 12 

j ei {c, z) = Kilh(c, z), (282) 

where ki is a proportionality factor determined by the normalization. 
It is desired that as cz — * °° , 

J e*(c, z) — > ^==- cos (cz — ^ 7r^ . (283) 

Thus it follows that 

ki = i l+1 x/ 2, (284) 

and 

Jei(c, z) = a/I E' i l ^D l n J n (cz). (285) 

IK 2 n=0.1 

Since U\ — K X V\ , one has 

j ei (c, z) = kjKx Se?{c, z). (286) 

12 The notation J2ej n , Re | 4) , Ro\ l \ Ro employed in an earlier paper (Proc. 
Nat. Acad. Sci., 6, 51, 1935) for the even and odd radial functions of the first and 
second kinds has been replaced by Jei , Nei , Jo i , Noi to eliminate a superscript 
and to emphasize a relationship with the corresponding solutions of Bessel’s 
equation. 
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Upon replacing a n by D n in (152), we obtain 

/Sej^Cc, z) = \Z2 tt \[ 6) Jei(c, z), 

where 

XS*> = *£V-'§, (ieve n) 

7T n“0 Do 

xS ” 1 = - Z' m-‘ ^ (iodd). 

C7T n“l Dl 


(287) 


(288) 


5.3 Even F unc tions of the Second Kind. The radial functions of 
the second kind are obtained by replacing J n (cz) in (285) by the Bessel 
function of the second kind JV„(cz). 

Ne t (c, z) = kiU & 

= a/ 1 ±'i l - n D l n N n (c,z). (289) 

2 n=0,l 

As cz — > 00 , 

iWc, z) -* ^ 7 ^ sin ^cz - ^ ^ • ( 29 °) 

The convergence of (289) is limited to the region | z | > 1. To ob- 
tain expressions valid in the neighborhood of z = 1 or z = 0 one must 
replace Lh by its equivalent expansion (198) in terms of hyperbolic 
functions. We shall define 

Se?\c,z) = Ua/Kx 



2 /«; 

xlJCj 

> 

/SeS 1} (c, z) H- E' [Dn sinh n rfi + \pD n cosh 

n=0,l > 

[ (291) 

where 






Dn = f +1 0 > °)» 

(292) 


ln- t + C - 
4 

n =«2 nDo 




i; .-{ D ; + -.)+i]j 

, (293) 


EVA 

n=0 
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for l even, and 

•n — 1 t \ 

Z Dn 

(n 2 — l)Di 

£ i -n( D : + 4(”-- 1 ) + J (294) 

r.“l 

for £ odd. C and a(x) are defined as in (195). 

It is obvious that 

Se i 2) (c, 2 ) = xS e) iVe/(c, 2 ) . (295) 

The Wronskian of Jei{c, cosh \p) and Nei(c , cosh \]/) is unity. 

Je, A Ne, - Ne, A 'Je l = 1. (296) 

# d\f/ 

Through combinations of these functions one may define Mathieu 
functions of the third and fourth kinds. 


f = + c + 2±' 

Ki 4 2 »=3 


Hei l) (c, z) = Jei(c, 2 ) + iiVe^c, z), (297) 

He { i\c , 2 ) = Jez(c, z) - ziVej(c, 2 ). (298) 

The notation corresponds to that of the Hankel functions. Similar 
functions can be constructed from Se ; (1) (c, z) and /Sel 2) (c, z). When 
cz — > 00 , 

1 . / 2Z+1 \ 

He l l) (c,z)-+^-e V' 4 (299) 

1 . / . 21+1 \ 

He?\o, z)->A=- ? y ez+ — V . (300) 

’ V C2 


5.4 Odd Functions of the First Kind. These are constructed from 
Ur- and Vr-i unctions subject to the condition a n = — (— l)”a_ n . Let 

F n = -iW^x a n . (301) 

Then according to (179) and (180) 

V 3 — —i ^ Fn sinh n\p, (302) 

n- 1.2 
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U _ _ K f\ fy ni~ n F n J ni.cz) • (303) 

We take 7, to be the standard V - or angle function of the first kind 
and define: 

oo 

iSoS l) (c, cos <p) = 22 f’nsinn^, 

n— 1,2 

normafized at z = 1, or tp = 0, such that 

fj-SoS^Cc, cos <p)1 = 1. 

\jLcp jip»o 

From this it follows that the F\ must be such that 


(304) 


(305) 


E = L 

n-1,2 


(306) 


For the standard U- or radial function of the first kind we choose 



JOi(c, Z ) = KlUzifi, z ). 

(307) 

The factor « is fixed by the condition that when oz -> ®, 



1 / 21 + 1 \ 

4 /' 

(308) 

whence 

= — i l nc/2, 

(309) 



(310) 

Since TJz 

= KiVz , one has 

JOl{c, z) = KZ-&3 5o? > (C, 2 ). 

(311) 

Upon replacing a* by F n in (165) and (169), we obtain 
iSoi^Cc, z) = \/2ir Az o) Joi{c, z), 

(312) 

where 

(1 even) 

XC 2 n-2 r2 

(iodd) - 

XC n-1 *1 

(313) 
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5.5 Odd Functions of the Second Kind. The radial function of the 
second kind is defined by 

Noi(c, z) = Kill 7 


? E 'nf-FMa), 

2 71-1,2 

(314) 

: -4 oo ; 


( 21 + 1 \ 

(“ - 4 y- 

(315) 


■ Not(c ’ 

The convergence of (314) is rapid only when cz is large. In the neighbor- 
hood of z = ±1 one must replace this expansion by its equivalent series 
(213) in terms of hyperbolic functions. Again by definition, 

SoS 2> (c, 2) = WKt 

2 K - 1 (316) 


where 


- , ^ so\ x) (c, z) - i jry cosh + F n’A sinh n $ f » 

7 r \Kz »-o,i 




(317) 


K's _ , c , 1 , 4 y> 

K,- U 4 +C + i + k (n*-4)ft 




(318) 


J^ni'F* 


n—2 


for J even, and 

As , o ^ 1 I 2 V' m ‘ + - 
K, = ln 4 + C “ 2 + h (n* - l)Fi 


fy r +1 F„ 


when l is odd. It is again apparent that 

iSof(c,z) = \/ 2 ^'k ( i o) Noi(c ) z). 


(320) 
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The Wronskian of the functions Joi(c, cosh \p) and Noi{c, cosh \{/) is 
unity. 

JoSno,- Not±Jo, = l. (321) 

d\p ay 

Odd Mathieu functions of the third and fourth kinds corresponding, 
to the Hankel functions are constructed in the usual manner. 


When cz — > 00 


Ho[ l \c,z) = 

Joiic, z) + iNoi(c, z ), 

(322) 

HoS 2> (c, i) = 

Jo t (c, z) — iNoiic, z). 

(323) 

*) - 

1 t(c* r r) 

V cz 

(324) 

Hof{c, z) - 

i ■ ( 21+1 \ 

l -*( c * — r ) 

-> — — e v 7 - 

V cz 

(325) 


Professor P. M. Morse and the authors are indebted to the Kennelly 
Fund of the Massachusetts Institute of Technology for funds to sub- 
sidize the computation of numerical values for the functions defined m 
this paper. In the case of the elliptic cylinder functions the coefficients 
D' n and F'„ , and the separation constants b, have been calculated for 
values of c over the range 0 to 4.5 at intervals of 0.2 and for I - 1, 2 
3 4 The separation constants for both prolate and oblate «P he ™ d 
functions have been computed for values of c over the range 0 to 5.0 at 
intervals of 0.2, and for m and l = 0, 1, 2, 3; m + l £ 3 The coeffi- 
cients dL likewise have been calculated for the same values. These 
tables eventually will be published. Until then they may b< i ob tamed 
in mimeographed form by application to the Department of Physics at 

M The'authors wish to acknowledge their indebtedness to Dr. R. Albagli 
Hutner who has verified all the formulas employed in this paper and 
who has supervised the computation of expansion “efficients and 
separation constants. Dr. Hutner contributed also the method de- 
scribed in Sec. 2.13 for obtaining the coefficients a n m the case of Math 
functions of the second kind. 


Massachusetts Institute of Technology. 
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' This introduction to the useful. 

•oritulas for *^**£o“ CO ntained in the preceding 
Chey are based o Stratton. For the saKe of 

article by J. Chum slzed . namely, that Chu 

parity one fac ^ coefficients of the 

Stratton, when discussing whlch occur 

spheroidal functions, deal with the* , £ d to the 

» « *s r:: T< - « is—— 

tables are the ones, » 


o * ^ ~ ^ 

a n = 1 d 


d* (n+2m)± 
n nl 


t?t.t.tpTTC CYLIH DER-SfflCTIfl^ 

” “ “ “ v 


2o2Tdlnat£5. s 
X = | cos cp cosh if 
y = | sin <p sinh if 

N< 

04 <p4 # 

04 y < 00 

d = distance between foci 

Wave ftguati nn yields, 

& * ( b* - 
dtp 2 

+ (c 2 cosh 2 ty-b£ )R =0 

, -p = s or R, each of the above 

If z - cos<p = cosher, an 

equations transforms into 

('l-z 2 ')F"-zF , +(bn-c z )F - 0 

( • „in the paper of Chu and Stratton, 

which is the form given in the pap 
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The Even Angular and Radi al_Functian& for a given value of b x (c), 
lo Angular Function 

o® , 

• ZL n! 


where 


Se^ (c.coscp) = , l D n cos ( n(p ) 

oa 

HV = 1 

n=0,l U n 

s e„ (c,l) = 1 


Primed summations are over even values of n if i is even* 
over odd values of n if 1 is odd- 

r o, & l 


P 2ii (o, 

/ Se n (c*cos<p) Se, (c,coscp)d(p = < 

J 0 1 K 

\ = *L 2(b o )2 + Sa (D “ )2 ] 


N = •*£.' (D*) 2 
^ n=l 


k-JL 

1 = even 
JL = odd 


2. Radial Functions. 

oo J 


Je (c fZ )=4 H 'os(c^) 

x- n-u,i 

QQ ? 

Ne (c, 2 )=4 H sin (c 

n— 0*1 


- 2m±l v 
A ' 


See Jahnke and Emde for a definition of N n (x), 

For z=cosh\jr equal or almost equal to one* the solution 
jfe^ may be evaluated by means of the following formulas; 
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J L =s even : 

Nej (o,coshf) 


/afEfi II’ D 1 ’ slnh nV + (L^E^cosh n* 

-i */ «\E t £o “ i i n= ° 


where 


n=0 


i 11 

d oM I H*£ T" 


V 


E, 


OO | 

”5 & ‘*i 

L D n +D n| 

A 

CT ( §-l) 
n <L 1 ) 



s.. ». (*) .. «* of « im '“" C,maer r “" 10 “ 

for evaluation of C, c(n) and D n 


f. = odd : 


3r{% Sl' S D ”“’ h n * 


Ne, 


where 


oc* 1 

.=n 


■n+l^X 


E = 4^ ni“ T *D. 
n=l 


n 


v 


d iH 4n( * l - 2 &x 


(nlzil - 


n+3 n C , -r / n-ll 

DiZ-, 1 ■W Dir( Tj 

n=l < 


E 


X 


E 


1 


See Note (A) at end of section on Elliptic Cylinder Functions 
for evaluation of C, cr(n) and D n . 
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3. Joining Factors. 


Se 1 (c,z) 


\ 



n/2# ^Je^C 0 * 2 ) 


is 

* n= 0 





°n 

5 


t= even 

X= odd 


4. Wronskian and Special Values of the Radial Functions. 

Jen (c.coshy)— %e n (c,cosh^)- Ne fl (c^cosht) Jen (c, cosht) 
X cL\|r * ^ dt x 


At z=coshty kl: 


Je fl (c, cosht) + — 
i Xj^ 

— Je« (e,coshty) >0 

d 4r * 


Ne 


cosh t) ^-n/2« 


— Ne 0 (c,cosh>jr) > JTH X. 

d^ ^ ^ 


The quantity ^ is evaluated by means of the second forms 
given above for the Ne^ solutions. 
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m m ..... .jMmw * !l ' , “ v * 1 ” ot b i Co '- 

l o AngulaT Function. 

QO , ? 

So (c,co S <p)=XL F^sin(ncp) 
n=l,2 


where 


K =1 

n=l,2 

— S0.(0,C0SCf)— ^ 

dcp ^ Y 


1 


( C,COS ^ So k(c^ cos ( f) d( P- 


0 , 
N| > 


M l 
k= i 


V' ? ,r,l \2 

n; = * 21 < F n ) 

x- n=l,2 


2. Radial Functions. 


oQ/ i „ o a .1 / 2 m+l .. 

, T ‘o \ * T fez) ^ COS(CZ- It 

,_fi J&L-L+ , 1 nF n J n toz; cz^- 4 

9 z)-*l 2 z n=l,2 


Jo A (c,^>-v 2 


= g JELYl ’ i n - < nF^ n (oz)^-^sin(cz 


lo (9»)-J5 z n=1>2 


2m+l 

4 


«) 


Solutions of HO^V) z = cos t equal or almost 
equal to one are: 


57 


L = even 


2 A 


NOl (c s cosh 


- 2 - U.A 

E 'i n=0 



where 




See 

for 


Note (A) at end of section on 
evaluation of C, <r(n), and F n 


Elliptic Cylinder 
* . 


Functions 


J?.= odd : 


cir 

(c j, cosh 10 = -j=Z 


5+1 


r !i£r F ^cosh n'|r+(L I + S)£Vsinh r 

H n=1 


F?t 


where 


y v ’ ,n+l v a 

\ m h. 1 n 


r -.“I «p r ' nitl .n+l fZ,!' + pi rS: 

H H )+c - i] +z E . x S L n * 




i = 


e / 


°¥>‘ 


See Note (A) at end of section on Elliptic Cylinder Functions 


Jhf 

for evaluation of C*cr (n) , and F r0 


3, Joining Factors. 


i So (c,s)= J2Z y Jo t (c,z) 

v ■ s £; $ *- — 

x ac n=2 F^ 

v = a_£ i n - 1 4- x= odd 

jI ac n-1 F^ 

4 . wronskian and Special Values of the Radial Functions. 

Jo £ (c, cosh t) J- *>£<«, cosh t) - N0 £ (c, coshf)^Jo t (c, 


At z= cosh t ** • 

Jo^(c y cosher) — ^0 

■=«» *>— 

No£ (c , cosh t) — y 

— Ho (c, cosh y) ->J2* I*' 

dty X. ■*- 

The quantity |i< is evaluated by differentiating the 
second forms given above for the No £ solutions and then 

letting z=l. 


coshty) 
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Note (A) 


C = 0.577215665 


(Jahnke and Emde) 


a(n)= 1+ \ + J + ; i o-(o) = 0 


To 


obtain D and F^* (the superscript, Jl , is omitted 


for the sake of simplicity): 


let 


D 

^n D. 


or 


n-2 


_n 

n-2 


then 


1 d + °n.2> 

7 n =* 2 -f 

* - n - ? Y n + 2 


where 


6 = b. - — 

l- 2 


n,2 


, n^2 

• n=2 


Yn " « 2 d««,2> n 


2n + 


^n+2 


d +a n +2,2^ 


2(n+2)+Yn +4 (2 [n+ 4 >' 


v» 

T n 


n 


D A 


n 


^n-2 

V2 


or 


y» F» F f 9 
T n _ n _ n-2 


Y F 
Y n n 


n-2 
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Q. g evens 


D& = 0 , D£ = D 0 Y^ 

= 1 > D i = I D i Y i 



PROLATE SPHEROIDAL-F UNCTION S 


Coordinates 

X= f sinh\jr sin© cos cp 
y= | sinh sin© sin cp 
_z= ~ cosh t cos 9 

Lets t] = cos 9 > 

d= distance between foci. 


z 
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Wave Equation Yields : 


frf ( n 2 -l) 5T + [ A+c ^ 2 - Ofhr] 8 = 0 (K) 

it- < £ 2 -d ft- + 2 - -fzr ] R = 0 (L) 


Let 

Z= T) 

= € 

m 

and 

S or R = 

(z 2 -!) 2 w 

then 

(l-z 2 )w" - (2m+l) 

zw*+ j^-A~m(m+l)-c 2 z^ j w = 0 

where 

-A-m(m+l) = b in 

the paper of Chu and Stratton 

Angular Solution. 




s<;>,(=,o> ■ c t <»> ■ #•" *' 


•q=0 


d Cn i 1 


dTl 


n=o 


4 "i (O’ n > = P m + l (T l ) 

primed summations are over even values of n if £ is even, 
over odd values of n if l is odd. 


| f* ' (a£) 2 (n+2m)! ^ g =k 
n^O.l nt(2n+2m+l)’ 


The summation of the angular function is from n 

— » » 1 *• - 

schematic representations , using 

Stratton. 

S U> (c.ti ) = (l-n 2 ) 2 i d n' I n (T1 ) 


'ml 


n 






where 





(Spherical Bessel Function) 




cz-)- 150 


cos (cz- 


X+m 

2 


*) 


where 

*W< Cz) 



N n+m+ 


1 (cz) = 

2 



The summation of these radial solutions are from n= 0 or 1 
rn=~ , as can be seen from the following diagram: 

In the notation of Chu and Stratton 




S’ 


-00 


(or-2m+l, for n=odd) n- 


finite 


+ oo 


W|H 


When the above formula for does not converge 

(namely, z=l or close to it), or converges too slowly (namely, 
m= large, such as 3), the following formulas can be used: 


^-= evens 



4 


(n+2m) l 
n! 




n 


dV. 

n^m+n v 


Jc = odd: 




-8c 


m-2 




l-2m 




^4C»< 


In these cases the summation extends from -oo to + o® « 


m 


N m 


^m+n W = T “ («) 


p is a quantity which approaches 0, as n approaches integer. 

(See Chu and Stratton) 



r 


Q=fin ite x p 


finite 


lx 

+ oo 


-2 m 0 

(or -2m+l,for n=odd) n 
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OBLATE SPHEROIMLJMCI1MS 



(M) 


(H) 


IYir wave Wnuatlon_ yigj^S s 

<n 2 -» *, * p" 2 -^l s ' ° 

' , +ic equation (U) becomes of the 

If c is replaced by ^ hence wi il have the same 

— f °~ “ r ° 0 :“ “ » ;P.> 1= » - 1 ' 

solutions, provided c kv 
in the solutions . 

If . 1. r=pl.~* W «'• -f ' “ 

« - *• is 

the solutions. 
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For example , a solution to equation (M) is; 


-tx3 t 



where the coefficients fjJ are, of course, different in 
value from the coefficients d^ of the prolate case, but 
are obtained in the same way. 


All the prolate equations for , B^j£, and 

apply for the oblate case provided that; 


1. c (prolate) be replaced by -ic (oblate), 

2. cosG =T\ remains the same, 

3. the coefficient d£ be replaced by f£, 

4. z=£ = cosher for the prolate case be replaced by 

z= = i sinht for the oblate case, 

5„ the positive square root is to be chosen in all 
cases . 


SPECIAL VALUES FOR THE OBLATE SPHEROI DAL FUNCTIONS 


At C= 0 : 

s S ( - ic ’ o) 


± 1 + m 2 m ^ 1 +2m-l_ ^ j 


(4) I 


=0 


„ , ji c m 4(2m) I 


=0 


/= 
i = 
1 = 

1 = 


even 


odd 


even 


odd 
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=0 


i 1 - 1 ^ c m ^ X ff (2n»l) I 


Ji = even 
1= odd 


n=l 


f(-i+2S=i-)! 

X = even 

2m+3 c m-2 (Aj£^~) I y 

R i 2 L ( - iC,i0) " 

Itl fi i n+1 ( Jif=au.) I + i £ 1 ’ f ^ n (2af=i) 1 (¥) 5 

2 n=l n J 2 n^l 

(|)1 

^ ' n-l £i ^Sli- t= Odd 

* (_i) »£. i 1 

. x 

5 =0 JS (2) * ( m “ 2^ lf -2m § 0 n nl 


ft R S ( ' ic,iC } 


5 b 1 ^ 


:g° f » f i & 

+ (_l) m+1 2 H i n ~f S 

v n=2m+2 V 2 ; 


l- 


even 


r c %)^ 2 


i ■ 1 1 
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_ rt-p r> and A (or c and B) das 
For certain values of c and A t numer _ 

been found, both theoretically and as a result of num 

ical calculations, that one or more of the negat 

coefficients for the range However, 

-2-l<»*l f0r B °“' TZm be shown that each 
_i -i ,. uc h cases investigated, it can De sno 

1 11 a ,,(1) s (2) (see Chu and Stratton), R^Z 

of the solutions S mj _ '- see m 

or rOO behave properly. The essential points of th 

argument are: (1) » is always an integer, (2) we wi 
ways observe the rule that p— >0 before c— >c k , where 
0 1 = a quantity which approaches zero as n approaches 
whole 8 number (see Chu and Stratton), and c, is the value 
of c for which one or more negative coefficien s 

Infinite. 


MOTES 

In addition to the formulas to be found in the 
first part of the article by Chu and Stratton, the 
following will be found useful: 

^ m - (j-ffl-l'fj. for |n|> 2m 

V-n P 

(to be used for evaluating part of the negative half 
of the summation in the second solution for R^)- 

(2) For m = positive integer, and n unrestricted: 

. nnm+n 


m 

n f.i ) m ? m+n (m+n) !. (n+2 m ^ ,( 2 

<+n (z)= 1 (2n+2m+lJ» z n+2m+l 


X 


F(a±2S±2, n±|S±l. m+n+ 1. ij) 
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Some auxiliary tables available, 
u Bessel Functions 

** Jn(l) ^ "Tables of Functions" 

a. Jahnke-Emde, T 

. tiTheory of Bessel Functions 

to. Watson, "Theory 

. of the British Association for 

e. Report of the 

the Advancement of Science, 

** NnCX) _ „ e n Tables of Functions" 

a, Jabnlee-Emde, T 

^ Watson, "Theory of Bessel Functxons 

Report of the British Association for 

the Advancement of Science, 19 H 

9 „ 


N. 


n 


a g 

x n 


lx J n+1 (X) and J. n -i W 

Jahnhe-Emde, "Tables of Functions" 
h_ Report of the British Associatxon 
the Advancement of Science, 1925 

. of th e British Association for 

Sx Report of the „ l6> 

the Advancement of Science, 

1922 ^ \ W 


j 1 w = £ Si 

n+2 


(*) = 


n 


JL 


C B W 


-n^ 


d . W.P.A. Tables to be published soon. 
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■^sbesi 


Bo, Legendre Functions 

p m (COS0) and its derivative 
a! Jahnke-Ende, "Tables of Functions" 
b, Tallquist, Acta Societatis Sciential 

Fennicae, Vol. 33, 1908, Ho. 9 
W.P.A. Tables to be published soon 

2, p» ( X ); ^(i*) “ d thelr deri ™ tives ’ f ° r x> 

^ W.P.A. Tables to be published soon 
^ Q * w> ^(ix) and their derivatives for *>1 

L W.P.A. Tables to be published soon. 


R 0 Albagli Hutner. 


7 ^ 




TABLES 



in the tables of the spheroidal coefficients, the 

— — s: rrr;:L“ . 

cant figure. Thus: _ 

.(3) 24692 = .00024692 

If the sign is the same throughout the whole length 

of a zzx .ua - ***>*- - the r p of the c 

otherwise « is given in front of each number. 

Although five significant figures have been published 
ta the tables of coefficients, the fifth 

considered doubtful, ^^.^^/^roidal Functions 

the row c-1.4 in . „ the calculations 

.,3, ,.o. «r ■”* 

— » - — 

A few of the sets of separation constants, k fnd B mf 

T~rzriz 1 ::.^ - ^ •— 

m=2 , J0 • 

in some of the tables of coefficients with ne ^ ative 

subscripts, the coefficients become infinite for a certain 
subscript , falls between two con- 

value of C, namely, c . When lnflni ties has 

secutive tabulated values of c, 

been inserted in the table. 
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